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CHAPTER i 
INTRODUCTION 
1.1 Overview 


Modern bridge designs often include decks whose spans are at a skew 
angle to their supports for economic or aesthetic considerations. The 
behavior of these structures is largely dependent upon a variety of 
factors such as the angle of skew, direction of reinforcement, aspect 
ratio, orthotropic stiffeners and types of loading. The effects of these 
and other factors on the mechanics of deformation may be interdependent, 
thus requiring a study including variation of parameters. This report 
will examine the effects of skew angle on major principal moments of free 
Spans under simulated truck loading. Orthogonal and skewed reinforcement 
will also be discussed. 

The finite element approach has been selected as the method of 
analysis. A mesh size study is undertaken in Chapter 2 to examine the 
effects of singularities caused by concentrated loads. GTSTRUDL [5] 
software is used for the stiffness analysis on a VAX 11-/50 digital 
computer. A service load simulator program is written in FORTRAN-77 
(Appendix A) and is run semi-interactively with GISFRUDL to incrementally 
adjust the truck position. The parameter studies are discussed in 
Chapter 3. Here the major principal moment at center span is plotted 
against the location of the truck for various angles of deck skew. The 
moment reductions which are found for increased angles of skew do not 


necessarily allow for a reduction in reinforcing steel. The angles at 





which the major principal moments intersect the steel plays an important 
role in the efficiency of resisting flexure. This relationship is 


discussed in Chapter 4. 
1.2 The Nature of Plate Bending 


An orthogonal isotropic plate subjected to transverse loading 


deflects according to the differential equation [18] 





In the case of a simple span such as a bridge deck of Figure 1.1(a) the 


solution to (1) must also satisfy the boundary conditions 


w= 0 (Zero Deflection) (2) 
— Along Simple Supports 
2 2 
e+ v ay = (Zero Moment ) (So) 
aX dy J 
and 
3M aw 
a + v a = 0 (Zero Moment) (4) 
OX dy 
Along Free Edges 
(aw 4 33 - 
==: (2 - v) ——s) = 0 (Zero Shear) (5) 
OX dX Oy 3 
where: w = displacement 
x,y = cartesian coordinates 
v = poisson’s ratio 
q = loading function 
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These equations can be solved for many loading conditions by the use of 
Fourier series. For skewed plate analysis, the introduction of an 
oblique coordinate system, such as that shown in Figure 1.1(b) is 
required for the analysis. After transformation of the Laplacian, 


equation (1) becomes [17] 





D ae a yy . a a ae rh 
mosh oat paaeie 2 sinad) aa + Asind Sy + —- + wt mea 
(6) 
where: x = x + y tand 
y =y sec¢ 
@ = deck angle of skew 


p(x, .y) = transformed loading function 
Solutions to (6) are difficult even for simple loading cases and are 
therefore not well suited to extensive parameter studies, although Kennedy 
[10] has been successful with series solutions aaa computer. Analytical 


methods have also been applied by Krettner [12] and Lardy [13]. Energy 


methods have been used by Guzman and Luisoni [6]. 


1.3 Numerical Methods 
Several researchers have been successful in analysis of skew slabs by 
the finite difference [8], [£9] and finite element methods [1], [15]. The 
growing popularity of numerical methods for skew slab analysis may be 
attributed to both the inadequacy of analytical solutions and to the 
advances in today's computer technology. The latter, particularly, has 
allowed for extensive use of numerical methods in applied mechanics. 


Several package programs such as STRUDL [14] are available to the public 


and have become a widely accepted method for analysis, design and research. 





The Finite Element Method was selected for this study using GISTRUDL 
[5] software on a VAX 11-750 computer. The finite element approach breaks 
up or "discretizes’" the structure to be analyzed into a network of consti- 
tuent elements. In plate bending each of these elements is usually 
allowed three types of displacements or "degrees of freedom" at each 
corner node (two orthogonal rotations and one transverse displacement). A 
mathematical function is assumed to describe the displacement variation 
between the nodes. Then, the stiffness properties of the individual 
elements can be developed in matrix form. If compatability between the 
adjacent elements is satisfied then the solution obtained should 
"converge’' as the number of elements is increased. Here compatibility 
refers to the fact that the pieces must fit together and that all 
adjoining elements at similar nodes must have corresponding degrees of 
freedom. .Strictly, compatibility is not completely satisfied for the 
Bending Plate Parallelogram (BPP) element which is used in this study. 
Zienkiewicz [20] shows it is not possible for a simple polynomial expres- 
Sion to ensure full compatibility when only one ‘neonacenane and two 
rotations are prescribed at the nodes. However, experience [20] with the 
BPP element shows that it "converges" to a good engineering approximation 
in most practical cases. 

Although many different types of elements may be used in discretizing 
a structure, the procedure is fundamentally the same. The material 
properties and boundary conditions are first defined for the problem. The 
element mesh is then selected and the structure stiffness matrix is formed 
from the known element stiffness matrix. Matrix algegra can then be used 
to solve the equation relating displacement, moment, strain and stress. 


Several texts (such as £3] and [20]) are available which detail the finite 


element process and give many diverse applications. 





CHAPTER 2 
Mesh oles Tou 
Zl elnibCOGuet LON 


The first step in analyzing a structure by the finite element 
method is selecting a mesh size and pattern. The mesh is a very 
important feature of the study and must be selected with care. In this 
Study the following aspects of mesh size, type and pattern were 
considered: 

1. Accuracy of results 

2. Economy of computer time 

3%3* Ease of comparison of results at varied angles of skew 

4. Assurance of conforming to the criterion of non-distorted 

elements. 

GTSTRUDL provides for a wide range of elements suitable for this 
Study. The element type chosen is the Bending Plate Parallelogram (BPP) 
which is very well suited to skew plate analysis. The BPP uses a fourth 
order transverse displacement expansion and uses three degrees of 
freedom (one displacement and two rotations) at each corner node. At 90 
degrees the BPP is equivalent to the more familiar Bending Plate 
Rectangle (BPR). The element is not considered distorted unless it is 


skewed to an interior acute angle of less than 30 degrees. 


Consequently, it can be used for deck skew angles of up to 60 degrees. 





This study will model the deck as a free span between simple 
Supports. The thickness is held constant at 18 inches and the moments X 
and Y are released at the joints along the supported edges. Before 
proceeding further it must be pointed out that there are two ways to 
geometrically "skew" a plate from a rectangle to a parallelogram, both 
of which are used extensively in the literature. The first way is to 
keep all edges the same length and allow the supports to move closer 
together as the angle of skew is increased. This is illustrated in 
Figure 2.1(b). The other widely used convention is to maintain a 
constant distance between the supports and allowing the free edge to 
"stretch" as a skew angle is increased. Figure 2.1(c) illustrates this 
convention. This study will adopt the first convention for the simple 
reason that a thin strip taken parallel to the inclined side will 
approximate the span of an equivalent simple beam, the length of which 
will be held constant as the deck skew angle 1s varied. Moreover, the 
reinforcement is generally laid parallel to the deck edges as well and 
holding these lengths constant may be a more realistic approach. Thus 
for all the examples in which the results are compared between the 
skewed deck and an orthogonal deck, the geometries shown in Figure 


2.1(a) and 2.1(b) are implied. 
2.2 Determination of Element Size 


Before skew angle can be investigated, the proper element size must 
be determined by a convergence study of the solution of rectangular deck 
of dimensions similar to those used in practice. Hereafter, the term 
"span" will refer to the inclined distance (skewed dimension) between 


supports. Also the ratio of width to span will be abbreviated WSR where 
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the width is the distance along the supported edge. The WSR parameter 
will be varied throughout the study in much the same way as the aspect 
ratio 1S in orthogonal plate analysis. The span length chosen for the 
mesh size study was 30 feet and the supported edge length 45 feet 

(WSR = 1.50). Three mesh sizes were investigated based on these 
dimensions. 

The first mesh was called "RIGHT" (Figure 2.2) and consisted of 24 
elements (4 in the span direction and 6 along the supports) each 90 
inches square. The idea of the study is to start with this coarse mesh 
and compare the results to increasingly more refined meshes. However, 
Since the comparisons are made based upon a 10 kip concentrated central 
load, the finer meshes will not be approaching a finite value (because, 
as is well known from thin plate theory, a point load produces infinite 
stress}.-- Therefore it is anticipated that the true solution will lie 
somewhere between the first course mesh run and one of the finer 
modifications. A more accurate finite element analysis can be made 
using the load area equal to the tire contact (imprint) area, and 
comparisons can be made with the "point load" mesh runs to determine 
which will give the best approximation. The second mesh was titled 
"RTFINE" (Figure 2.3) and was given 8 span elements each 45" and 10 
transverse elements each 54". For the third run these 80 elements were 
bisected bilaterally to give a 320 element mesh called "RTEXFINE" 
(Figure 2.4). As expected the values for the central span moment under 
the load increased in the two successive finer meshes (see Figure 
2.5). The values went from 3.25 k-in/in (RIGHT) to 3.86 k-in/in 


(RTFINE) to 4.48 k-in/in (RTEXFINE). 
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2.2.1 Convergence 

Tne problem now is to determine which (if any) of these results are 
close to the moment found under a more realistic tire load. The double 
wheel contact area of a water tanker truck was studied in Reference 4. 
The dimensions are approximated here to a rectangle of 10" x 24". 
Interior element loads (i.e., concentrated loads within the boundaries 
of an element) can be approximated by placing the force upon imaginary 
Stringers parallel to the element edges and using the resultant reac- 
tions at the nodal points as a new collection of equivalent loads. This 
technique will be used later in the program SKEW LOADER as a method for 
finding statically equivalent nodal loads for any truck position on a 
skewed deck. This technique, however, is not applicable at this point 
as the objective here is to avoid the singularities caused by concen- 
trated -loads altogther. Therefore a more exact approach 1s required 
before we can obtain a solid basis for comparison in choosing the 
appropriate mesh size. 

Since GTSTRUDL will accept uniformly Peenibeted loads only over an 
entire element (i.e., no partial element loading), an extremely fine 
localized mesh has to be developed to surround this tire load. Extra 
care must be taken here to ensure that compatibility conditions are 
satisfied between element and that no elements are distorted. 

Distortion is defined in GTSTRUDL as having aspect ratios greater than 
2.0 or acute angles less than 30 degrees. Rectangular elements cannot 
be used in this localized mesh as compatibility cannot be maintained 
without highly distorting element aspect ratios. Therefore, a 
triangular pattern was carefully assembled for the localized mesh 


(Figure 2.6) using GTSTRUDL element Dending Plate Hybrid Triangle 
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(BPHT). Due to the hiyh degree of correlation between the data obtained 
from three different meshes (RIGHT, RTFINE and RTEXFINE) for areas not 
in the vicinity of the load, the RTIFINE mesh was chosen to house the 
localized triangular mesh. In other words, Figure 2.6 was inserted into 
the hatched portion of Figure 2.7. The central portion of the localized 
triangular mesh contains four rectangles each 5" x 12" which are loaded 
with a uniform pressure of 0.04167 ksi. This modified system is 
therefore statically equivalent to a central 10 kip concentrated load on 
the plate. This mesh, called UTRTFINE, gave a central span moment of 
4.107 k-in/in which falls slightly above the moment given by the RTFINE 
mesh (see Figure 2.8). Since the moment of the RTFINE mesh with a 
concentrated central 10 kip load is only 6.2% lower than this "exact" 
value, it is considered that the RTFINE mesh will be accurate enough for 
the purpose of investigating the variation in moment with deck skew 
angle. The cost in computer time for using the RTEXFINE mesh would be 
increased exponentially as the number of elements are quadrupled, and 
would as equally overestimate the moment as the RTF INE mesh would under- 


estimate it. 
2.3 Effect of Aspect Ratio 


Before the RTFINE mesh could be selected as the proper mesh size 
for the study, one more check had to be made. If the results showed 
substantial changes in accuracy when the aspect ratio of the plate was 
varied, then this would indicate that the mesh size used in the study 
would have to vary as the WSR 1s changed. A Similar run was therefore 
made to compare the cutput of the RIFINE mesh with a more exact value 


where the aspect ratic has been reduced from 1.50 to 0.60. The 
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respective meshes were renamed RTFINEO60 and UTWSRO6O (Figure 2.9). The 
mesh of Figure 2.6 was renumbered to be housed in the mesh of Figure 2.9 
as shown and comparisons were again made between the central span 
moments given by the two outputs. The results showed no substantial 
change in behavior at the new aspect ratio. The UTWSRO60 (more exact) 
gave a central span moment of 5.69 k-in/in while the RIFINEO60 gave 5.44 
k-in/in, or 4.3% less (see Figure 2.10). The RTIFINE mesh is therefore 
considered a satisfactory approximation to the true deck behavior and 


will be used as the basis for investigating the effect of skew. 
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CHAPTER 3 
PARAMETER STUDIES 
a eeelntiroducuron 


Now that the mesh RTFINE has been selected as a guide for the 
element size, the effect of skew angle can be studied. In order to 
obtain a realistic evaluation of the variation in maximum moments, the 
loading used must correspond to the dimensions and axle weights of an 
actual truck. Furthermore, a consistent convention must be used to 
establish the truck position for varying angles of deck skew. The 
principal moment at key points of the slab will be evaluated by GTSTRUDL 
finite element analysis for each position of the truck as it moves 
incrementally across the span. Then the corresponding results for 
different angles of skew can be compared when plotted on the same 
graph. Since the WSR may affect the moment variation as the skew angle 
is changed, it too will be treated as a parameter. 

Figure 3.1 shows the models to be studied. The principal moments 
will be computed for WSR's of 1.50, 1.00 and 0.75 while the skew angle 
is changed from 0 to 40 degrees. For the case of WSR = 1.00 (rhombus), 
skew of 20 degrees will also be investigated. In all cases the truck 
will move along the inclined centerline as measured from support to 
Support. In the case of the rhombus, the effects of edge loading will 
also be examined as the truck is moved across the span inset at 48" from 
the left free edge. In all cases the deck is simply supported at the 
top and bottom edges and free along the left and right edges. This will 


be the conventional] orientation of the deck throughout this report. 
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The truck chosen for the investigation was the 70 kip FDOT-Su4 
(Figure 3.2). The span dimension of the deck was selected for 
convenience to be 32 feet in order that 8 span elements of 48" (Figure 
3.3) may be used as in the RTFINE deck. The RTFINE actually used 8 span 
elements of 45", though the element study showed the difference should 
be of negligible order. The transverse dimensions were varied from 48 
to 32 to 24 feet to allow the WSR to change from 1.50 to 1.00 to 0.75, 
respectively. Transverse element dimension was chosen as 48" for 
convenience. Again, these parallelogram element dimensions remain 
constant as the skew angle of the deck is varied (Figure 3.4(a)). Note 
the reference position of the skew angle in Figure 3.4(a) as the 


complement 1S sometimes used in the literature. 


3.2 Procedure 


tee hae? 


For the span chosen for the study, the maximum centerline moment in 
an orthogonal deck would occur under axle #3 (wheels #5 and #6 in Figure 
3.2). Wheel #5 was therefore chosen as the reference from which the 
position of the truck will be measured. The truck will always be 
positioned as if it were moving parallel to the free edges of the 
deck. Since the program SKEW LOADER (see Appendix A) allows for input 
in terms of skewed coordinates, the position of wheel #5 will be given 
in such a coordinate system. The X input will be aS measured from the 
left edge of the deck and the Y input will be the span (inclined) 
distance from the base support. For center loading, the position of 
wheel #5 will be input such that the center of axle #3 1s over the 
transverse centerline of the deck (see Figure 3.4(b) for an example). 


The truck will be "moved" alony the span direction while maintaining 
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this centered position. In this way the advantage of symmetry can be 
maintained on the right deck and axisymmetric loading can be maintained 
on the skewed decks. 

The following discussion applies to the analysis procedure for load 
cases A and B (see Figure 3.1) where the WSR is 1.50 and the skew angle 
is 0 and 40 degrees, respectively. The remaining cases are handled in 
an analogous manner except for the edge loading cases which will be 
discussed separately. The first position of the truck is located at the 
transverse center and 96" (i.e., 2 elements) up from the base support. 

A finite element analysis is then run with the truck in this position 
and the principal moments at joint 59 (deck center) are obtained. The 
truck is then moved up 12" (i.e., one quarter element) in the span 
direction to a new position 108" from the base support, finite element 
analysis is again run with the truck in the new position and the 
principal moments at joint 59 again noted. This procedure is repeated 
at increments of 12" until the truck is 96" (i.e., 2 elements) away from 
the far support. The largest principal moment at joint 59 1s plotted 
against the truck position in Figure 3.5 where the letter label of the 
graph corresponds to the respective case in Figure 3.1. 

The above procedure is repeated for each of the remaining cases and 
the results are plotted for varied angles of skew at each WSR to obtain 
Figures 3.5, 3.7 and 3.9. In this way the effects of skew can be 
inspected separately for each WSR. Now for each of these graphs the 
transverse variation in moment must be studied on neighboring nodes in 
order to be sure that the moment at the transverse center is actually at 
or near the maximum in the deck. This is due to the fact that as deck 


skew angle is increased, thc element edges are moved closer to the 
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wheels of the truck which may allow for greater localized moments at 
joints off-center. Thus each of the Figures 3.5, 3.7 and 3.9 showing 
moment variation with truck position is followed by a Figure 3.6, 3.8, 
and 3.10, respectively, showing the section moment variation across 
midspan in the vicinity of the peak moment. This section diagram 
examines the moment variation at midspan across two elements to the left 
and right of the peak and therefore plots for a total of 5 joints. For 
example, the joint 59 moment-position graph of case A (WSR = 1.50, 
Skew = 0.0) peaks when wheel #5 is positioned 192" from the base 
Support. Therefore to investigate the transverse variation in moment 
for this position of the truck, Figure 3.6 plots the moments at joints 
41, 50, 59, 68 and 77. This graph is labeled "A-192-RT" giving the load 
case, truck position and deck skew angle, respectively. 

For the case cf edge loading a similar procedure is followed, 
though edge loading 1s considered only for the cases with a WSR of 
1.00. The truck is again placed 96" from the base support but now 48" 
(i.e., one element) in from the left free edge. A finite element 
analysis iS run with the truck in this position and the moment.at joint 
5 is noted. The truck is then moved across the span in increments of 
12" with a finite element analysis run for each successive position. In 
each case the moment at joint 5 is noted, and as before the truck is 
Stopped at 96" from the far support. The moment at joint 5 (the edge 
node at midspan) is plotted against each position of the truck and is 
shown in Figure 3.11 for skew angles of 0, 20 and 40 degrees. Ina 
Similar procedure as for the center loading cases, the transverse moment 
variation is plotted in Figure 3.12 which examines the moments across 4 


elements starting from the left edge. Fur this parameter study, 
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including both the center and edge loading cases, a total of 170 finite 


element runs were made. 
Soo EResulEs 


In many of the skewed decks, the maximum moment does not occur at 
the same joint as in the corresponding orthogonal deck. However, it is 
seen that the increase in moment from the corresponding joint never 
exceeds 3%. Therefore the moment at the same joint will be used for 
comparison. It must be noted however, that the moments generally do not 
peak for the same position of the truck. For example, Figure 3.9 (WSR = 
0.75, center load) shows that the right deck attains maximum moment at 
joint 32 when the truck position (i.e., wheel #5) is at 192" from the 
base support. The 40 degree deck however peaks earlier when the truck 
position is at 144", This is the general trend in all center load cases 
Studied here. The right deck always peaks at 192" while the 40 degree 
deck peaks about 48" (one element) sooner. As might be expected, the 20 
degree deck peaks about 24" earlier. ; 

The fact that the moments peak earlier for tncreased angles of skew 
can be attributed to the effect of wheel #6 (also a critical wheet) 
reaching the deck center sooner. This allows the moment at the deck 
center to peak earlier and also accounts for a portion of the reduction 
in maximum principal moment from the right deck because both wheels 5 
and 6 are not lying along the span center at the same time. fhe right 
deck peaking at 192" is to be expected as this corresponds to the 
position where the critical wheels are at span center. For the case of 
the edge loading however, it is noted that the joint 5 moment peaks 


occur when the truck is positioned at 192" cegardtess of the angle of 
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Skew. Figure 3.12 shows that the principal span moment decreases from 
the left edge of the deck at low angles of skew, though it increases 
Slightly at joint 14 at 40 degrees. This would seem to indicate that 
the localized effects of the free edge are less pronounced at greater 
angles of skew. 

Table 3.1 summarizes the key data for the parameter study. It is 
interesting to note that the greater WSR's allow for the greater 
percentage in principal moment reduction for increased angles of skew. 
For example, it is seen that the WSR of 1.50 has a moment reduction of 
24.3% in moving from the right deck to 40 degree skew deck. This can be 
compared to only a 15.6% reduction for the WSR of 0.75. The rate of 
moment reduction 1s also non-linear as evidenced by the table entries of 
WSR = 1.00 (center loading). Here the moment reduction for the first 20 
degrees (i.e., 0 to 20) is only 5.4% while for the second 20 degrees 
(i.e., 20 to 40) is 13.5%. Therefore the rate of moment reduction is 
increased with an increase in skew angle. In the case of edge loading, 
the behavior is similar except that the percentages are cent 
increased. It is important to note that these decreases in maximum 
principal moment do not necessarily allow for commensurate reduction in 
reinforcing steel. The angle at which the reinforcement is laid and the 
direction of these principal moments play an important role in the 
amount of steel required. These effects will be discussed in greater 


detail in Chapter 4. 
3.4 Contour Plot Description 


Appendix C shows the contours of major principal moment and Z 


(transverse) displacement for skew angles of O and 40 degrees. These 





T’€ S3YNI4s WOU 
uV¥u JISVO SALVOIGNI [V] 








G4IGNLs 
LON 


AYVWWIAS SAIGNLS YALIWWYVd T° JTVL 





aL 0 


00°T 


Ost 


JNIGVO1 d4LNso 














42 


plots were obtained from GITSTRUDL graphics and plotted on a Tektronix 
print device from the scope environment.! The loading conditions 
Selected for comparison were a central 10 kip load and the peak SU-4 
central load. The WSR is 1.50 throughout. 

Figures C.1 and C.2 show the effect of skew angle on the major 
principal moments under a central 10 kip load. Note that the contour 
gradient remains essentially perpendicular to the supports. This same 
effect can be seen on Figures C.3 and C.4 (SU-4 peak loading case) where 
the path of "steepest descent" is nearly along the shortest line to the 
supports. The concentration of contours near the obtuse corner show 
that the major principal moment increases more quickly (along a line 
towards the center) than at the acute corner. Figures C.3 and C.4 can 
be compared to the transverse moment variation graphs for load cases A 
and B as shown in Figure 3.6. 

et: C.5 through C.8 illustrate the variation in deflected shape 
between skew angles of 0 and 40 degrees. Figures C.4 and C.5 refer to 
the central 10 kip load while C.7 and C.8 refer t6 the SU-4 peak load. 
It is seen that the behavior is similar for both loading cases. As 
would be expected physically, the displacement gradient is perpendicutar 
to the supports and the direction of "minimum descent" is along the span 


centerline. 


1 The scope environment is a characteristic operating domain in 
GTSTRUDL which allows for interactive graphics at the terminal CRI. 





CHAPTER 4 
REINFORCEMENT 
4.1 Introduction 


The safe and economic proportioning of reinforcement is critical in 
bridge deck design. Ideally, the re-bars should be laid orthogonal to 
the principal design moments. However since the principal angles vary 
from point to point throughout the deck for even a single loading case, 
the reinforcement cannot be placed ideally in a practical sense. 
Furthermore since design is often based on a series of different loading 
cases, the principal angle will often vary at the point as well. 
Therefore the concept of ideal reinforcement for a bridge deck is a 
trivial one. 

Certain general directions for the reinforcement are a yood deal 
more economical (in terms of required steel PEC TE onty) than 
others. Cope [2] has studied orthogonal reinforcement for skewed decks 
and compared the experimental results of placement parallel to the 
Supports versus parallel to the skewed edge. The results of this study 
will be discussed in more detail in section 4.3. Orthogonal 
reinforcement in skewed decks however, has severe limitations in 
practice. Since a large number of re-bars of different lengths are 
required for this design, a great deal of extra labor is required for 
cutting and placing the reinforcement. Therefore, though it is 
generally desireable to have reinforcing steel running perpendicular to 
the supports (to resist the span moment), such an arrangement ts usually 


not practical. 
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Discussions with FDOT engineers indicate that most reinforcement in 
skewed decks today 1S not orthogonal. The primary reinforcement is 
generally laid parallel to the deck edges and the transverse reinforce- 
ment is laid parallel to the supports. This is certainly the easiest 
from a construction point of view, but for largely skewed bridges it may 
require a considerable amount of extra reinforcement. Morely [16] has 
developed a design process for non-orthogonal reinforcement in skewed 
decks. Though the method may be tedious for practical design problems 
witout the aid of computer, it serves well to illustrate the effects of 
skewed reinforcement. This design will be discussed in Section 4.4 and 
some examples worked which are relevant to the parameter study of 
Chapter 3. 

Even though the STRUDL analyses are purely elastic, it is relevant 
" first. examine the failure mechanisms of skewed decks by the theory of 
yield lines before discussing reinforcement in detail. For a more 
complete discussion on yield line theory and applications, several texts 


are available such as [7], [19]. 
4.2 Failure Analysis by Yield Lines 


Analysis of plate capacity can be determined theoretically by the 
principle of virtual work and an assumption of a yield line pattern. 
The work done externally by the applied forces acting through the plate 
deflection is equated to the internal work done by the rotation of the 
moments acting on the yield lines. The collapse load can then be solved 
for in terms of the ultimate unit moment capacity of the plate. Often 
the main problem in the yield line approach is determining the correct 


(or critical) failure mechanism to be analyzed. Yield line patterns can 
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be assumed from some experience or from some generally established 
guidelines. Hughes [7] gives some such guidelines, although generally 
the critical pattern is not easily found even for some rather simple 
geometries. 

For example, Hughes L7] shows three possible failure mechanisms for 
a simply supported skewed span as shown in Figure 4.1. The deck is 
given a central point load W and the orthotropic reinforcement is 
oriented parallel to the deck edges as is generally used in practice 
(see Figure 4.1). The figure and equations have been modified to agree 
with the conventions used in this report thusfar. It is seen that the 
elliptical fan pattern represents the least allowable failure load. 
Since bridge decks are generally not subjected to a single central 
concentrated load, this case may be of purely academic interest. 
However, in deriving the critical load for the elliptical failure 
mechanism, Hughes outlines (and proves) a very useful principle for 
obtaining an affine isotropic right slab from an orthotropic skewed 
one. This procedure is discussed below. 

A skewed deck such as that shown in Figure 4.2(a) can be 
transformed into an equivalent orthogonal deck for the purposes of 
analysis. The equivalent right deck (or so called "affine deck") shown 
in Figure 4.2(b) is of course much easier to analyze and can be 
obtained by the following rules: 

(a) Deflections are identical at corresponding locations in both 

decks. 

(b) Given that mand um are the ultimate resisting moments in the 


reinforcing direction of the actual deck, then the affine deck 


has ultimate resisting isotropic moment m. 
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(c) Given that the Y coordinate axis is along the m reinforcement 
direction in the actual deck, then all distances measured in 
the Y direction are the same for both decks. 

(d) Given that the X coordinate axis is along the um reinforcement 
direction in the actual deck, then it is taken at a right 
angle to the Y axis for the affine deck. 

(e) In order to obtain a dimension in the X direction of the 
affine deck, divide the corresponding length on the actual 
deck by \fu. 

(f) Divide the loads on the actual deck by \/u cos to obtain the 
corresponding loads on the affine deck. 

Using the above procedure, an interesting result is obtained if the 
deck in Figure 4.2(a) is analyzed assuming a yield line occurs across 
midspan- (see Figure 4.3). This assumption for the failure mechanism may 
be valid for the multiple load cases arising on bridge spans. The load 
P is a sum of the concentrated wheel loads acting across midspan and 
should cause failure to occur in beam action. If m represents the unit 
moment resisting capacity of the affine deck in Figure 4.3(b), then by 


Statics 


(4) (a Pacg) = my) 


eee 
4b cosé 


Therefore the deck of Figure 4.3(a) is equivalent to an orthogonal deck 
of span length L and base of a length diminished by the cosine of the 


angle of skew. This illustrated in Figure 4.4. 
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4.3 Orthogonal Reinforcement in Skewed Decks 


An experimental investigation was conducted by Cope [2] on the 
effects of orientation of orthogonal reinforcement in 45 degree skewed 
decks. Although orthogonal reinforcement is generally not used in 
skewed bridge design, the study finds some interesting results which may 
lead to a better understanding of flexural behavior. The test procedure 
was fundamentally a comparison between two 45 deyree skew slabs, each 
with orthogonal reinforcement. Slab A had secondary reinforcement 
parallel to the supports. Slab B had primary reinforcement parallel to 
the skewed edges. The propagation of cracking and deflections in the 
two slabs were examined for increasing loads in various locations. 

The tests showed considerable behavioral differences in the two 
Slabs in the areas of deflections, cracking and failure mechanisms. 
Sagging cracks were first initialized on slab A running towards the free 
edges and parallel to the supports, and under heavier load hogging 
cracks appeared at the obtuse corners. Initial sagging cracks in stab B 
occured earlier (at about 80% of the load for A} as did the hogging 
cracks at the obtuse corners (at about 67% of the load for A). The 
cracking patterns were similar in direction and spacing, though somewhat 
Straighter and more continuous in slab B. The modes fo failure of the 
two decks were different as well. Slab A initially fatled in shear on 
the free edge at the obtuse corner. At about 40% higher load, stab B 
developed a top surface crack between the two obtuse corner supports 
which was wide enough to produce a discontinuity in slope. Both slabs 
were able to carry considerably more load after their initial 
failures. Ultimate failures were punching shear at the obtuse corner 


for slab A and excessive deflections for slab B. 
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The report concluded that the behavior of skewed decks is strongly 
influenced by the direction of the reinforcement. When the reinforcing 
is placed orthogonally and parallel to the supports, the slab is stiff 
and behaves well under service loads. However, heavy concentrated loads 
are not distributed well across the slab and a large reactive force 
develops at the obtuse corner. This led to local failure in the test. 
For slabs with orthogonal reinforcement parallel to the free edyes, a 
more flexible slab results which better distributes moments due to 
concentrated loads. Greater deflections can occur and ultimate failure 
load 1S increased, though again hogging cracks at the obtuse corner may 
limit servicibility. 

The results of this study are consistent with a finding by Kennedy 
[11] which experimentally investigated the stress near corners of simply 
Spanned.skewed plates. Here too it was found that stresses near the 
obtuse corners are Significant and increase with increasing angles of 
Skew. Furthermore, the stress at the obtuse corners may exceed the 
maximum stress at center span. Kennedy then ree that the obtuse 
corners of concrete saved decks be heavily reinforced top and bottom in 
directions parallel and perpendicular to the supports as well as 
parallel to the free edge. The accute corners should be heavily 
reinforced at the bottom and nominally reinforced at the top. The 
direction of reinforcement near the acute corners should be 


perpendicular to a line between them and parallel to the free edge. 
4.4 Design for Skewed Steel 


A design procedure was worked out by Morely [16] on proportioning 


skewed reinforcement which is laid parallel to the edges of the deck. 


ne 


ne 


| 
7 





3 


Since this is generally the orientation used in construction, this 
method could prove useful in design practice. Equations are developed 
based on the optimum proportionment of steel to resist a moment triad at 
a point in the deck. Therefore the method begins with the assumption 
that elastic analysis has been performed for a load case on the deck and 
that a set of three moments (two bending and one torsional) are known at 
the point considered. The design procedure is then carried out using a 
series of charts, tables and equations given in the reference. 
Morley 's design for skew reinforcement to resist a single moment 
triad consists of the following steps: 
(a) Orient the sign convention of the known moment triad in 
accordance with Figure 4.5 
(b) Compute Be wey an k = (ele using a rough 
SS igaeton for 19. (k=0 if no minimum steel is specified) 
(c) Enter the charts with the nearest values of $ and k and locate 
the appropriate region 
(d) Evaluate aj,, Aig? 42x» Fag» oe and Or, eon the appropriate 
equations in the tables 
(e) Compute 1, by an iteration equation of section equilibrium 
(f} . Compute Aj,, Aig? hae Ao? from aj, = ot aAny (x «> 6) 
where these variables and others used in the design examples are defined 


as follows: 


Coie specified minimum steel are per unit slab width 
Aix» ayy ‘area functions’ with dimensions of moment 

da» do depths for bottom and top steel, respectively 

k failure factor connecting aan and applied 


twisiny moment Myy 
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l lever arm of tensile steel in failure direction 

ae My» Myy bending and twisting moments per unit slab width 

Eee} principal moments es 

Ty, To total resolved steel forces in failure direction 
at bottom and top, respectively 

O.> BF optimum failure directions for sagging and 
hogging 

> skew angle of reinforcement 

Oy, yield stress of steel 

OO value of concrete stress block at failure 


A total of five examples will be worked here to examine the effects 
of skewed reinforcing steel. The first three examples show the extra 
reinforcement required for a 20 degree skew deck over a right deck when 
designing for the central span moment for a 100 kip load. Examples 4 
and 5 similarly illustrate for the peak moments produced by SU-4 
truck. It should be noted that this design procedure could prove very 
tedious for practical problems without the aid of 3 computer. This is 
because the moments at many points in the slab would have to be 
considered for various combinations and positions of vehicle loads. The 
results for each trial run are very sensitive to the sign and magnitude 
of the torsional moment at the point considered. These multiple moment 
triads must be considered and the minuimum steel retained at each 
successive iteration. Morley makes allowances for this in an extension 
of the procedure for multiple moment triads along with some guidelines 
for a rational approach to design. 

It should be understood that in the examptes which follow, the 


steel in the slab is not campletely designed for the given loading 
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condition, since only the moment at the geometric center is considered in 
each case. The examples are intended only to illustrate the effects of 
non-orthogonal reinforcement at a point for a few different load condi- 
tions. Morely emphasizes the effects of the sign of the torsional moment 
by working two similar examples of a 20 degree skew deck. The first has 
a point subjected to My = +50 k-in/in, M, = +35 k-in/in and Myy = +20 
k-in/in and the total required steel is found to be 3.96 in¢/ft. In the 


Second example the sign of M,,, is changed to -20 k-in/in and the total 


5 
amount of stee! becomes 8.20 in¢/ft, r roughly twice the amount. This 
increase occurs because the direction of the major principal moment is 


shifted, due to the sign change of M from within the acute angle 


xy? 

between the reinforcement to the obtuse angle. This reduces the unit 

effective reistance of the steel, and hence a greater amount 1s required. 
For the examples discussed below, reference is made to Figure 4.6 


and the following data is applicable to each: 


Deck thickness = 18.0 in 
E = 3324 ksi : 
fo = 3401 psi 
o. = 2891 psi 
WSR = 1.00 
oy, = 50 ksi 


In addition, the following equations from Morley [16] are common to all 


examples: 
_ 1 2 Ne 
li Ty = y— lai, cos OF yy Sin (0 - 6)] 
a 
re ee «die (12) dee (rar 
: a e 1, O 2 O. yy, 


iene ay, =o. 158 (x «> $) 





Se eS 
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Example 1. Right deck 100 kip central load 


Step (a) GISTRUDL gives the following moments in the appropriate sign 


convention 

= 28.49 k-in/in 
My = + 43.03 k-in/in 
Myy = 0.0 k-in/in 

M M 


oat et 
Mey xy 


(e) 


(f) 


Using Figure 2 [16] the appropriate region is A 
From Table 1 [16] 
doy = aay = 0 + No Top Steel Required 


ajy = Ys + Myy | = 28.49 


ayy = My + IMyy| = 43.03 
0. = -45 
From equation 1, 1, = +~ [28.49 cos(-45) + 43.03 sin“(-45)] 
a 
i = Some sal = 0 
4 
From equation II, 1, = 16 - 202-891 


or 


2 Z E 
1,¢ - 16 1, + 6.18 = 0 > 1, = 15.60 


d 

From equation III, 12(28.49) = 50(15.60)A,, 
Aly = 0.44 in/ft 

12(43.03) = 90(15.60)A; 


Ayy = 0.66 in@/ft 


yy 
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Example 2. 20 deg skew deck using moment triad of example 1 


Step (a) 


(b) 


(c) 
(d) 


aye C849 k = 
My = + 43.03 k-in/in 
Myy = 0.0 k-in/in 
M M 
a ees 
Ty TP | xy! 
Using Figure 3 [16] the appropriate region is A2 
From Table 1 [16] 
2x = aog = Q + No Top Steel Required 
: = 2 (1 +2 sin 6) * 
= My tan cepa UY tang - Mey 
aj, = 28.49 - 43.03(tan 20)* + 2 +2 sin 20) [4303 tan 20-0 
, ayy = "50255 
M 
y 1 
ay, = oe te (M:,Cttano - MO 
1$ cos“é cos $¢ ' y Xy 
a4 4 = eee entice {43.03 tan 20 - 0| 
(cos 20) 
diy = 65550 
2 ll 
0. = x (90 + 4) 
ig 0 
os oo 


From equation 1, Ty = + [50.85 (cos 20)° + 65.36 (sin 35)*~ 
a 


ia = 66.4/1.5 To = 0 


, 
| : 1 
From equation II, 1, = 16 - a75-a5F 
or 
2 2 
lao - 161, + 11.48 = 0 


I = Seco" 





(f) From equation III, 12(50.85) 


Aix 


12(65.36) 


a 


Comparison of this result with Example 


skew reinforcement minimized for the same moment triad. 
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i 


50(15.25) Aj, 
0.80 in2/ft 


I 


50(15.25)A) , 


1.03 in¢/ft 


1 illustrates the effects of 20 deg 


The total steel 


per unit width becomes (0.80 + 1.03) cos 20 = 1.72 in¢/ft or 56% more than 


Example l. 


Example 3. 


20 deg skew deck 100 kip central load 


Step (a) GTSTRUDL gives the following moments in the appropriate sign 


convention 
=a = 2496 K-10 AN 

My = + 40.77 k-in/in 

Myy = + 0.646 k-in/in 

M M 

Kee) Os y _ 40.77 _ 

(c) Using Figure 3 [16] the appropriate region is A2 
(d) From Table 1 [16] 

a2x = aay = 0 > No top steel required 

7 2 (1 + 2 sing) 2 
ayy = My - My tan“¢ +e IM tan 4 Ie) 

9 1 +2 sin 20 _ 
aly = 24.36 - 40.77 (tan 20) + (a5) [40.77 tan 20 0.646 } 
aiy = 44,39 

M 

_ z 

a OS IM tan ¢ ra 
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f ee + ye 140.77 tan 20 - 0.646| 
aig = 61.24 
8. = +5 (90 +4) 
on =e 
(e) From equation 1, 1, = ~~ [44.39 (cos 55)* + 61.24(sin 35)°] 
T) = 34.7/1, . i ee 
From equation II, 1, = 16 - DIDBITY 
or 
1,° = -161, + 6.0 = 0 
1, = 15.62" 
(f) From equation III, 12(44.39) = 50(15.62)A), 
Rie On6e inc/tt 
2 12(61.24) = 50(15.62)Ay 
Ary = 0-94 in@/ft 


Comparison of this result with Example 1 illustrates the additional 
Steel required at the geometric center for the 20 deg skew deck vs. the 
right deck, both under a central 100 kip load. The total steel per unit 
width becomes (0.68 + 0.94) cos 20 = 1.52 in*/ft or 38% more than 


Example 1. This increase is due to the major principal moment lying 


within the obtuse angle of the reinforcement as shown in Figure 4.7. 
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Right deck peak SU-4 loading 


Step (a) GTSTRUDL gives the following moments in the appropriate sign 


convention 


M, = + 5.07 k-in/in 


= 
1 
ae 


13.90 k-in/in 


—= 
I 
a) 


-0 k-in/in 


M M 


pi : 
i ay | 


(c) 


(d) From Table 1 [16] 


Ady = a2y = 0 > No 


a Aix = My + IMyy | = 5,07 


= 13.90 


aly My + IMyy| 
6 = 


: -45 


From equation I, Jy 


Using Figure 2 [16] the appropriate region is A 


Top Steel Required 


9.49/1,, Tp = 0 


q 


From equation II, 1, = 16 - 5757-g97y 


or 


is - 


1, = 15-90" 


(f) From equation III, 12(5.07) 
Alx 
12(13.90) 


Aly 


| 


= =]61, = 1.64 = 0 


50(15.90)A,, 
0.08 in@/ft 
50(15.90)Ayy 


0.21 in@/ft 
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Example 5. 20 deg skew deck peak SU-4 loading 


Step (a) GTSTRUDL gives the following moments in the appropriate sign 


convention 

M, = +4.66 k-in/in 
My =  +13.13 k=in/ io 
Myy = -0.36 k-in/in 


M 
m4 4.66 _ y « Weed su. 
a)” oss” Tht” 0.367 = “°° 


(c) Using Figure 3 [16] and noting that M 0 


Xy . 


M, > -4.66 , My > -13.13 


and the appropriate region is B 
(d) From Table 1 [16] with suffixes 1, 2, s, h interchanged 


a9x = agg = 0 + No Top Steel Required 


_ 9 (et esr) 
ajyx = ~My A: My tan=9$ Eos oe 


| 


Mt _M 
ene - Ry 


Me col 1aiaii(tan 20)° = (15,13 tanmz0) + 0-36] 





COs 
aly = 10.84 
eet) | 
a = + M tan o-M™M 
1o cos<$ cos $ ' y xy 
13.13 1 
Wiese? ° cos 20 |-13.13 tan 20 + 0.36| 
(cos 20) 
ary ey 
By) _ £60 
@ = ty (90 +) = 55 
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i 
i 


dq 
Ty = 16.0/1,, To = 0 


ye + 19.57 (sin 35)*] 


(e) From equation I, Ty [10.84 (cos 20 


; 
From equation II, 1, = 16 - PBST) 
or 


2 = = 
wy CU LON, + eo = 0 a = 15.82 


(f) From equation III, 12(10.84) 50(15.82)A,, 


Ay, = 0.16 iné/ft 
12(19.57) = 50(15.82)A, 4 
= inZ 
Ayg = 0-30 in?@/ft 


Comparison of this result with Example 4 illustrates the additional 
Steel required at the geometric center for the 20 deg skew deck vs. the 
right deck, both under peak SU-4 loading. The total steel per unit 
width becomes (0.16 + 0.30) cos 20 = 0.43 in¢/ft or 49% more than 
Example 4. This increase is again due to the major principal moment 


lying in the obtuse corner of the reinforcement as shown in Figure 4.8. 


» 





66 


SNIGVOT b-NS AVdd 
SJ IdWVX3 LNAWSDYOSNIZY 8° 3YNDIS 


p JidWVX3 YIAO GaYINdIY SI 
TA3LS WLOL JYOW %6p =“ SUWa-3y 
JHL N33ML39 JTONY ASNLEO FHL NI SLNSWOW TWdIONTYd HLIM 
S317 WdIONIYd YF9NV1 = S 31dWvXd (q) SSQIONIOD 1351S % AIdWvX3 (2) 


NOTLOSuIO 
uvd-3sy 





(UL/UL-y) O6°EL = We W 








CHAPTER 5 
CONCLUSIONS 
5.1 Summary 


The finite element method was used to examine the effects of skew 
angle on the major principal moments of simple spans. A mesh size of 
eight span elements in 30 feet was found to give results within about 6% 
of those found from a more intricate localized mesh surrounding a more 
realistic tire load. These results varied little with change in aspect 
ratio. This mesh size was therefore chosen as a guide for the parameter 
Studies which followed. 

Skew angles of 0, 20 and 40 degrees were investigated for various 
width to span ratios (WSR's) as shown in Figure 3.1. The major 
principal moments at key points in the deck were examined as an FDOT 
SU-4 type truck was moved across the span. The service load simulator 
program of Appendix A was used to calculate the equivalent nodal forces 
for the truck in any position on the skewed deck. The results showed 
that for center loading, peak moment reductions of up to 24% were found 
for 40 degree skew decks over orthogonal decks. The percent reductions 
were lower for smaller WSR's. In general, the major principal moments 
peaked earlier (i.e., at lesser advanced positions of the truck along 
the span) for increased angles of skew, This can be attributed to the 
fact that the wheels of the axles do not reach span center 


Simultaneously. 


6/7 
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The direction of reinforcement plays an important role in the 
flexural behavior of the deck [2]. Orthogonal reinforcement laid 
parallel to the supports provides for a stiffer slab which behaves wel] 
under service loads. Orthogonal reinforcement laid paratlel to the free 
edges provides a more flexible slab. However, skewed reinforcement 1s 
the preferred method from the construction point of view. The affine 
deck method discussed in section 4.2 seems to be the most practical for 
use in the design office. Morely's [16] method for minimum steel at a 
point would seem to have limitations in practice, although it serves 
well to illustrate the considerable effect that the direction of the 


major principal moment has on the quantity of skewed steel required. 
522 Furtnes study 


More extensive parameter studies could be used with different truck 
types and combination loading. Program SKEW LOADER (Appendix A) could 
be easily modified for this application. In addition, the effects of 
girders, deck thickness and material properties Zaita also be studied. 
The reinforcement design method developed by Morely [16] is well suited 
to computer programming. This could be an area of further study along 
with complete design examples and comparison with the affine deck 
method. 

The stress concentration effect at the obtuse corners could be 
investigated by the finite element method. This could be of particular 
importance as experiments have shown that this is often the area of 
initial failure. An extremely fine local mesh could be assembled (such 


as along the line of Figure 2.6) in the area of the obtuse corner and 


the moments investigated with variation in deck angle of skew. 





APPENDIX A 
PROGRAM SKEW LOADER 





APPENDIX A 


PROGRAM SKEW LOADER 
A.l Introduction 


Since the parameter study calls for moving a truck incrementally 
across a skewed deck and evaluating the moments for each truck position, 
a method had to be developed for loading the deck. GTSTRUDL allows for 
concentrated loads to be placed only at the nodal points of the mesh. 
Thus the weight of each truck wheel within an element must be broken 
into its statically equivalent nodal forces before it can be input to 
GTSTRUDL for stiffness analysis. To do this by hand would be extremely 
tedious as 4 reactions would have to be calculated for each of 8 wheels 
for 170 different load cases. Furthermore, the positions of the wheels 
in relation to the mesh nodes are geometrically complex making large 
numbers of calculations prohibitive. This type of problem is therefore 
well suited to computer programming. 

Program SKEW LOADER was developed to output the equivalent nodal 
forces for a truck positioned on a skewed deck. The program will 
generate any mesh size on a deck of any dimensions at any angle of skew. 
The data for the truck FDOT SU-4 is stored internally in the program. 
Although any type of truck may be used if the truck width, axle spacing 
and axle weights are known. The user also has the option of specifying 
the truck position in rectangular or skewed coordinates. The latter is 
recommended for most applications and in generat greatly simplifies the 
input procedure. The position of any wheel may be input to establish 


the truck position, which is always assumed to lie parallel to the left 
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and right edges of tne deck. The program uses the stringer method with 
each loaded element to solve for the equivalent nodal forces. 

A set command and regeneration prograns were also written to 
expedite the study. These programs were used to operate on the first 
GTSTRUDL program by rewriting the joint load input portion based on the 
output of SKEW LOADER. In this way the truck could be repositioned 
automatically and sent back to GTSTRUDL quasi-interactively. This 
provided for an orderly and systematic approach to data acquisition. 
This program set was written in command language on a VAX 11-750 
computer. The program SKEW LOADER is written in FORTRAN-77. 

The rest of Appendix A gives a flowchart, program description and 
input guide for SKEW LOADER. Copies of the program are also included 


along with two examples and sample outputs. 








A.2  FLOWCHART /2 


| 


INITL 


LOCATE 


GMESH 
















TRUCK 

TYPE FDOT 

SU-4 
? 






READ IN 
TRUCK DATA 





READ IN PHI, LCN, WNUM, XSHIFT, YSHIFT, 
XINPUT, YINPUT, KEY, NXSTEP, NYSTEP, 
XINC, YSKINC 


GENERATE SKEWED MESH STORE 
NODAL COORDINATES IN XCOR, YCOR 





ae - 
ro bee 





FIND COORDINATES OF LEFT FRONT WHEEL 
BASED ON COORDINATES OF USERS 
CHOOSEN INPUT WHEEL : 










COMPUTE COORDINATES OR ALL WHEELS 
BASE ON POSITION ON LEFT FRONT WHEEL 












ARE 
ALL WHEELS 
ON THE 

DECK? 





PRINT WARNING FOR EACH 
WHEEL OFF THE DECK 
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~ ~| pO FOR EACH WHEEL 


FIND WHICH ELEMENT THE WHEEL FALLS IN 
BY SEARCHING FIRST IN THE SPAN THEN IN 
THE TRANSVERSE DIRECTION 


4 
| 
| 







SEARCH 


COMPUTE THE FOUR NODAL POINT REACTIONS 
FOR THE ELEMENT BY USING 
IMAGINARY STRINGERS 


on ee ee ces eee me ee ee ee ee ee ee ee ee ee ee oe 


mse: KEEP RUNNING TOTAL OF NODAL 
WHEEL | REACTIONS IN MATRIX REACT (I, J) 








PRINT OUT NON-ZERO NODAL POINT 
REACTION WHICH GIVES THE 
EQUIVALENT DECK LOADING 







RESULT 
PERFORM STATIC CHECK 
TOTAL NODAL LOAD + WT OF OFF WHEELS = 
TRUCK WT 
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A.3 Program Description 


A.3.1 Main 

The program SKEW LOADER is structured as a series of subroutine 
calls from the main. The initial matrix sizes are dimensioned 
arbitrarily in the opening statement allowing for a maximum mesh size of 
30 nodes each way and a maximum truck size of 18 wheels. These numbers 
are only set for convenience so that variably dimensioned subroutines 
can be used later. Therefore the maximums above can be increased i 
by adjusting the sizes in the opening dimension statement (no other 
program changes are required). The input is read from LOADER.DAT and 
the output goes to QUIPUT.LIS. Only the key subroutines will be 
discussed below as the functions of the others are self evident. The 
subroutine descriptions below contain a complete listing of the 
(aniebles used within. Variables which are not listed are defined 


within the subroutine in terms of the ones described here. 


Aeoee Subroutine INITL 
This is where most of the user supplied data is read in, stored and 


echo printed. The following variables are used here: 


VARIABLE DESCRIPTION 

AXLSP Axle spacing matrix 

AXLWT Axle weight matrix 

DDX Deck dimension along X-axis 

DDYSK Deck dimension along Y (skew) -axis 
ERR Error switch 


KEY Mesh print request switch 


= ae 
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VAR TABLE BESCRIPTION 

LCN Load case number 

NA Number of axles 

NAXLSP Number of axle spaces 

NW Number of wheels 

NXCROS Number of crossings (mesh lines) X direction 
NXSTEP Number of steps (elements) X direction 
PHI Angle of deck skew 

TRTYPE Truck type (SU-4 or custom) 

WIDTH Width of truck 

WJX X coordinate of user's input wheel 

WJY Y coordinate of user's input wheel 
WNUM Input wheel number 

WT - Wheel weight number 

XINC Increment (element size) X-axis 

XSKINC Increment (element size) Y (skew)}-axis 
XSHIFT Shift coordinate en X-axis 

Monit | Shift coordinate switch Y-axis 


A.3.3 Subroutine GMESH 

The mesh is generated here based on the user tnputs specified in 
INITL. The nodes are generated first in the span direction and then along 
the X-axis. The coordinates of each node are stored in column vectors XCOR 


and YCOR. The following variables are used here: 
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VARIABLE DESCRIPTION 

ERR Error switch 

KEY Mesh print request switch 

NXCROS Number of crossings (mesh lines) X direction 
NXSTEP Number of steps (elements) X direction 
NYCROS Number of crossings (mesh lines) Y direction 
NYSTEP Number of steps (elements) Y direction 

PH] Angle of deck skew 

XCOR Nodal X coordinate storage matrix 

XINC Increment (element size) X-axis 

YCOR Nodal Y coordinate storage matrix 

YSKINC Increment (element size) Y (skew) -axis 


Petes routine LOCATE 

Now that the position of one of the truck wheels is known (user's 
input) the rest can be determined assuming that the truck is lying along 
the span direction (1.e., parallel to the deck's free edges). From 
trigonometry the position of the left front wheel can be found from the 
known position of the input wheel. Then a standard procedure can be 
followed to locate all wheels relative to this position by geometry. It 
may often occur that one or more of the truck wheels are lying outside the 
boundary of the deck. In this case a warning is printed to show the user 
that the total weight of the truck is not on the deck. This must be taken 


into account in the statics check at the end of the program. 
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The following variables are used in LOCATE: 


VARIABLE DESCRIPTION 

AXLSP Axle spacing matrix 

DDX Deck dimension along X-axis 

DDYSK Deck dimension along Y (skew)-axis 

ERR Error swtich 

NA Number of axis 

NAXLSP Number of axle spaces 

NW Number of wheels 

NXCROS Number of crossing (mesh lines) X direction 
NYCROS Number of crossing (mesh lines) Y direction 
OFFWT Matrix containing weights of wheels off deck 
PHI Angle of deck skew 

PSD Perpendicular support distance 

REACT Matrix of nodal reactions 

USERX X input (carried) 

USERY ~ Y input (carried) 

WHEELX X coordinate of wheel under consideration 
WHEELY Y coordinate of wheel under consideration 
WIDTH Truck width 

WJ X X coordinate of user's input wheel 

WJY Y coordinate of user’s input wheel 

WT Matrix of wheel weights 

W1X X coordinate of wheel #1 

W1Y Y coordinate of wheel #1 

XSHIFT Shift coordinate switch X-axis 

Wonk | Shift coordiante switch Y-axis 


3.5 Subroutine SEARCH 
This subroutine is the heart of the program and certainly the most 
complex. The function here is to find the location of each wheel with 
respect to its surrounding nodes, then compute the nodal point reactions by 


the stringer method. For a right angle deck the search is a relatively 
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simple procedure as the nodal coordinates can be easily compared to the 
coordiantes of the wheel. The skewed deck however, provides for 
complications in locating the position of the wheel relative to its 
neighboring nodes as a series of offsets from the edge of the skewed 
element must be computed and compared. The idea is to first search the 
Span direction for the first node with a Y coordinate larger than the Y 
coordinate of the wheel. Calling the Y coordinate of this node and the one 
immediately below it Dl and D2, respectively, the Tower "track" D2 can then 
be searched. Now the first node along this track with an X coordinate 
larger than the X coordinate of the wheel is noted. The X coordinate of 
this node is called CHECK1. Now the X offset distance of the wheel from 
the element edge is examined to see on which side of the line it falls (see 
Figure A.1). If XOFF1l is greater than WXOFF1 the wheel is in the element 
B: the-right. If it is less than WXOFF1 then it is in the element on the 
left. If the two distances are equal then it is on the element edge. A 
Similar procedure is used to locate the neighboring nodes above the wheel 
(j.e., On upper track Dl). Provisions are made a cases in which the 
wheel falls directly on a node or on an element edge. 

Now that the coordinates of the four neighboring nodes are known, the 
distances to the wheel can be easily computed. An immaginary grid of 
Stringers is laid on the element as follows: 

1. One is positioned under the wheel and laid horizontally just 

reaching the element edges 

2. Then two more are laid under the first along the skewed edges 
The reactions of the four corners from the weight of the wheel are then 
computed by statics and a running total of the results is kept in a storage 


matrix for the mesh called REACT. The following variables are used in 


subroutine SEARCH: 
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VARTABLE 


I,L,N 

K 

LCOUNT 
NA,NB NE NW 
NNE ,NNW NSE ,NSH 
NW 

NXCROS 
NXSTEP 
NYCROS 
NYSTEP 

PHI 
RA,RB,RL, RR 
REACT 

RNNE ,RNNW ,RNSE ,RNSW 
WHEELX 
WHEELY 

WT 

XCOR 

XINC 

YCOR 

YSKINC 


Rese) Subroltine RESULT 
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DESCRIPTION 


Node counters 

Wheel counter 

Element counter 

Nodes above, below, east and west 

Node northeast, northwest, etc. 

Number of wheels 

Number of crossings (mesh lines) X direction 
Number of steps (elements) X direction 
Number of crossings (mesh lines) Y direction 
Number of steps (elements) Y direction 

Angle of deck skew 

Reaction above, below, left, right 

Matrix of nodal reactions 

Reactions node northeast, etc. 

X coordinate of wheel under consideration 

Y coordinate of wheel under consideration 
Matrix of wheel weights 

Nodal coordinate storage matrix 

Increment (element size) X-axis 

Nodal coordinate storage matrix 


Increment (element size) Y (skew)-axis 


The non-zero nodal reactions obtained from subroutine SEARCH are 


printed out here and summed to check statics. The weight of the truck is 


checked against the sum of the total load on the deck and the weight of the 


"off wheels". 
the truck. 
VARTABLE 


NCOUNT 
NXCROS 
NYCROS 


The user can verify this result to be the original weight of 


The following variables are used here: 


DESCRIPTION 


Node counter 
Number of crossing (mesh lines) X direction 
Number of crossing (mesh lines) Y direction 





OFFWT 
REACT 
TOTAL 
TRKWT 
USERX 
USERY 
WJX 
WY 
WNUM 
XSHIFT 
YSHIFT 


A.4.1 


Data deck 
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Total weight of wheels off the deck 
Matrix of nodal reactions 

Running sum of reactions for statics check 
Back summed weight of truck 

X input (carried) 

Y input (carried) 

X coordinate of user's input wheel 
Y coordinate of user's input wheel 
Input wheel number (user's choice) 
Shift coordinate switch X-axis 
Shift coordinate switch Y-axis 


A.4 Input Guide 


The data file must be titled LOADER.DAT and be accessible from the 


main. ~The following data must be input in the order shown beginning in 


column #1. 


CARD 
] 
2 


FORMAT 


Free 


El 


BELO | 
1I2 
eZ 
i 


The term "card" refers to one line of data in the file. 


DESCRIPTION 
Alpha-numeric title of output 
Enter truck type as follows 

(a) 1 is FDOT SU-4 

(b) 9 is custom* 
Angle of deck skew (degrees) 
Load case number (user's option) 
Input wheel number 
XSHIFT switch 

(a) 1 is on 


(b) Ois off 








10 


11 
eZ 


13 
14 
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111 YSHIFT switch 


(a) 1 is on 


im} 0 1s off 
LE HZ X input for wheel given on card #5 above as 
follows: 


(a) if XSHIFT is on, give the X 
distance of the wheel from the 
left skewed edge of the deck 

(b) if the XSHIFT is off, give the X 
coordinate of the wheel 

late? <2 Y input for wheel given on card #5 above as 
fol lows: 

(a) If YSHIFT is on, give the distance 
to the wheel as measured along the 
inclined span edge for the base 
support 


(b) if YSHIFT is off, give the Y 


coordinate of that wheel 


Ol Mesh print request switch 
(a) lis yes 
(b) O is no 
eZ Number of elements in the X direction 
WZ Number of elements along the Y (skewed) 
direction 
Meteor c Element size X direction 
LEGe Element size Y (skewed) direction 
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*Note: If the custom option is chosen, then the following data must be 
entered on the next data cards (i.e., after card #2): 
(a) Number of axles [format 111] 
(b) Truck width [format 1F6.2] 
(c) Spacing for each successive axle as measured from the 
previous axle (one per data card) [format 1F6.2] 


(d) Axle weights (one per data card) [format 1F6.2] 


A.4.2 Examples 

The above proceure is illustrated by two examples. The first is an 
input to SKEW LOADER for the SU-4 case, and the second is for the custom 
option. Figures A.2 through A.5 given the input data and resulting truck 
position diagrams for the two cases. Comments are provided on the sample 


mpues fOr clarity. 


A.5 Extension of the Program 
The program can be easily extended to include a number of standard 
truck types beyond the SU-4. The easiest way would be to write more 
Subroutines in the form of TRTYPE1 to include data for any trucks. Then of 
course modify the switches in subroutine INITL. The program could also be 
easily extended to include more than one truck on the deck at a time. 
Finally, a more sophisticated mech generator could be included to allow for 


elements of different dimensions. This feature would be particularly 


useful in the analysis decks with stiffening girders. 
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¥ LOADER EXAMPLE 1 (SU-4) 


= Ae Rete Ate SY ED ED SD aD CE SD GD ane See cates eoes GS cutee atte OE ery SE SO eee Grey SOE ete em OE ete ete ete mete eet es GD i ae te SE cts Guy GED Ge GSD GD Oe 


cK TYPE 
-E OF SKEW (DEG) 
) CASE NUMBER 

JT WHEEL NUMBER 

YPUT (SHIFT ON) 

WPUT (SHIFT ON) 

4 PRINT REQUEST 


FDOT SU-4 


40.0 
1 
a 


192. 900 
240.00 


NO 


3ER OF ELEMENTS 
MENT SIZE 
{ DIMENSION 


TWARNING#### “WHEEL i 


FWARNING#### WHEEL # 2 IS OFF THE 


tt Orr 


48. 00 
384. 00 
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THE DECK TO THE NORTH ### 2WARNINGH#S# 


HHHHHHRHKHHKHKHKHKHHHKHHEHHHRHRHEHHHEREKE 
HHHEHKHKHKHKHKHKHKHHKKAKRHHKHRKREREEDHEE 


+ % 
# RESULTS OF LATEST ANALYSIS * 
+ + 


HHEHHKHHKKHKHHHKKHHHKAKAHKEHERREREE 
HHEKKHHHHKHHHHHKHKEHRERHHHHHEHHKEE 


SEL # WEIGHT 


9.39 
9.39 
Fa30 
9.39 
7.39 
9.39 


Ww N WD OI we OW 


RT O 


IN ELEMENT # 


oo 1 


vO 
ABOVE 29 
Do 


RT QO 


Cc 


re 


F 28 


4 


AFFECTING NODES 


ONLY 
52 63 
ONLY 
61 62 
OMLY 
50 61 


43 
71 
42 
A9: 
40 
&9 


44 
ya 


ae 
AL 


vo 


X COORD 


See 
63.81 
S727 
Goa 2 
659.937 

Leo. 07 


DECK TO THE NORTH ##*#WARNINGH iE 


Y COGRD 


222.15 
271.00 
183.85 
232./0 
145.99 
194.40 





SILTANT NODAL LOADS 


rH LOAD (KIPS) 
C 0.39 
1 8.96 
; 9.35 
c 8.96 
‘ 0.39 
( Gee 
8.36 
§. 36 
we QO 
: 0.45 
0.50 
0. 60 
Cae 


TICS CHECK 
AL LOAD ON DECK = 56.10 KIPS 
CK WEIGHT = 70.00 KIPS 


EKNOTE XXX 
CE THE XSHIFT IS ON X INPUT OF 192.00 GIVES X COORD 


HNOTE##%# 
CE THE YSHIFT IS ON Y INPUT OF 240.00 GIVES Y COORD 


37.73 


183.85 
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| LOADER EXARIPLE 2 (CUSTOM) 
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as + enaD Cnnee CY CED OAR CD CRNED NS SS ne ee i mie APE SEE Opt OES MED OOD OUPET Ore wees Cues ne tes ee one cami St ED OEESD OUDSY ONDSD Ott embaD Sti SEES apts HEPES orane 


OS A ee ee LT ES DS Cmte ORLED FEET OTLEY Orne SFOEP CUED GEIED CRLED CbSD ODED GEMED ceime cutee fmts emis amie Sie eetey ne? OEESY emime caine cmims amin @mbay OUED cwums quis ouums wuPeS amine enue? mes 


JK TYPE 
NUMBER GF AXLES 
TRUCK AXLE WIDTH 
AXLE SPACE # i 
AXLE SPACE # 2 
AXLE WEIGHT # 1 
AXLE WEIGHT # 2 
AXLE WEIGHT # 3 

-E OF SKEW (DEG) 

) CASE NUMBER 

IT WHEEL NUMBER 

PUT (SHIFT ON) 

IPUT (SHIFT ON) 

1 PRINT REQUEST 


CUSTO?P! 
2 
72.00 
168.00 
180. 009 
3. OO 
32. OO 
32.00 
30.0 
1 
3 
72.00 
192.00 
NO 


ied ep TS See empey Stet Ete au GREED CGD CHEST CEE) aus GENGD GEbGD GEbGD GHbED GENE GEbED SEE? GENE cue au tr creme cums GED cums orm om emms See cue cee cee cee cee cue CSSD freee Ore cee cee fee ee cee cee & 8 ee Se 


= ome omeee tie comer caper See que ae GS ees cess G6 cue eee Stet ember cscs SE GSS GSP CEP quem qemm fee SS em ce ee ee ee en eS ee me 


IER OF ELEMENTS 
IENT SIZE 
. DIMENSION 


13 ™ 
24.00 
350. 00 


WARNING#### WHEEL # 2 IS OFF THE DECK TO THE NORTH ### #WARNINGHSE# 


HHEHHKHEHERHHHHHHHKHHHKHHREEHEEEHE 
HHEHEEHHHKRHHHKHKRKRKHKHHHKRHHHKEEHKHERE 


52 t+ 
# RESULTS OF LATEST ANALYSIS # 
+ # 


HEHHKKKHHKEKHKHHRHEHKKEHKRHERREEHE 
HHEHHHKHHKHKRHHKHHHKEKHHKHHHHHREHREHEE 


EL + WEIGHT PMSECEMNENT 


4.00 ABOVE 3O 


AFFECTING NODES X CCORD Y COORD 


ONLY 48 ~108. 00 kok at 7, 





3 14.00 ABOVE 23 ONLY 41 

; 14.00 70 74 75 90 91 
: 16. 00 RT OF 14 ONLY 33 34 
4 14.00 63 67 68 83 984 


JILTANT NODAL LOADS 
j= + LOAD (CKIPS) 


: 8.00 
8. 00 
14.00 
4.00 
8.49 
i teye 
2.96 
8.09 
3.80 
ienlis 
1.33 
3.62 


B1CS CHECK 
AL. LOAD ON DECK = 68.00 KIPS 
2K WEIGHT = 72.00 KIPS 


FNOTEX## 
2E THE XSHIFT IS ON X INPUT OF 72.00 GIVES X COORD 


ENOTEX#H# 
ve THE YSHIFT IS ON Y INPUT OF 192.00 GIVES Y COORD 


—-24. 00 
38.359 
56. 00 

128. 35 


-24. 00 


1466. 28 


9] 


166.28 
202. 28 
[eo 
16.39 
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CHHEEHHEHEHHHHEEEE EHH E SERA REHEARSE RHEE 


KHEEKKEEH HEH HK HHHKHEEKE Ht HE 
HEEKHEREKEH HH HE KERKEREHR HE HH 
HR KEE HH EH tH et 
KHEEREKHE HKHRHEK RHEHHH a HX 
HHEHEEHHHH aat HH KHREKER HH * HE 
HE Ht HH He tH KH t- 
HHHHHHHE HT 5 ea FRHEEKKEKS He HH BH HEH 
HHRERHHHHH Ee EE KHEHEEHHHK Ee HH 
HH KEERHEKH EE HHEEHHG KEEHEKHE HEHHEE 
HH KEEEPHEH KEE KHHLELHH RHEHREREKE HEHHRE 
HE Kt te HH EE HE tie HE tt +H 
HH ea + HHRHHHHHH Kt KH HHEEREKE HHHEHRE 
HH EH HH KRKREEKEKEH EE HH HHEEKHE HMEHEHH 
eH HH Hk H+ Et HH RE we KH tc 
HHEEEHHH HKHEEKEHH ate HH KHEEERES HHEEHHHKEE EK 1035 
HHEEERRH HHEEHRHH HH HE KHHHEHHH HHRERKREH ooo. Et 


HEKHHHRHHKHHKEKKHKHEKKRRKHRRHHHKRRKHHEKRKRERKREHEHEEHRRRHEEHRRERSE SRE HRHEHBEERHER 


THIS PROGRAM CALCULATES THE SET OF EQUIVALENT NODAL LOADS FOR 
ANY TYPE OF TRUCK MOVING ACROSS A SKEWED (PARALLELOGRAM) DECK 
OF ANY DIMENSIONS WITH ANY MESH SIZE. THE MESH IS GENERATED 
INTERNALLY AND CAN BE QUTPUT GPTIONALLY. THE USER MUST INPUT 
THE FOLLOWING: 
(1) ’ALPHA’. GIVE ON THE FIRST DATA CARD THE ALPHA-NUMERIC 
TITLE 10 BE PRINTED ON THE OUTPUT LISTING. £77 CHARACTER 
(2) “TRTYPE’. ENTER HERE THE INTEGER CORRESPONDING TO THE 
TRUCK TYPE 
FOR WHICH THE ANALYSIS IS TO BE PERFORMED. THE OPTIONS 
ARE LISTED BELOW 
“™ (A) 1 IS FDOT TYPE SU-4 
(B) 9 IS CUSTOM 


F THE CUSTOM OPTION IS CHOSEN, THEN THE FOLLOWING 
ae MUST BE ENTERED ON THE NEXT DATA CAMRDS : 
) 
) 


WS AE AK IK KK A KR KK ae eK OK OK ae He FE OK RR a OK OK IR A Oe OK OK 


NIOOINAAAAAHANGANAAIAIAAIAIAIANANANAAANAAANANRANNND 


| 
| 
| 
| 
{ 
| 
| 
{ 
{ 
{ 
{ 
| 
{ 
{ 
| 
{ 
{ 
{ 
{ 
{ 
{ 
{ 
{ 
{ 
‘ 
{ 
{ 
{ 
{ 
{ 
{ 
{ 
| 
{ 
| 
{ 
{ 
i 
l 
i 
i 
{ 
{ 
{ 
{ 
| 
{ 
| 
{ 
{ 
| 
i 
i 
| 
| 
| 
t 
{ 
{ 
| 


ee ok ea ok ok oe oR oe ok ok ok 


HHENOTERHH 


NUMBER OF AXLES. CINTEGER UP TO 9] 
TRUCK WIDTH (COAXIAL WHEEL SPACING). CREALI 
SPACING FOR EACH SUCCESSIVE AXLE AS MEASURED 
FROM THE PRECEEDING AXLE (ONE PER DATA CARD). 
MOTE THAT THE # OF AXLE SPACES IS ONE LESS THAN 
THE NUMBER OF AXLES. CREALI 
(D) AXLE WEIGHTS (ONE PER DATA CARD IN KIPS). CREAL) 
ALSO NOTE THAT IF THE CUSTOM OPTION IS IN EFFECT THEN THE 
NUMBER OF DATA CARDS WILL VARY. 
(3) ’PHI’. THE ANGLE OF SKEW (DEGREES). CREALJ 
(4) SERECER LOAD CASE NUMBER (ONLY FOR USER’S REFERENCE) 
N 
(5S) “WNUM’. HERE THE USER GIVES AN INTEGER INDICATING WHICH 
ONE OF WHEELS HE WILL LOCATE. THE PROGRAM WILL AUTO- 
MATICALLY LOCATE THE REST (SEE SKETCH FOR RELATIVE | 
POSITIONS AND NUMBERING OF THE WHEELS). NOTE THAT THIS 
CARD COMPLETELY ESTABLISHES THE POSITION OF THE TRUCK 
ASSUMING THAT IT IS TRAVELLING PARALLEL TO THE CURB 
(SKEWED EDGE OF THE DECK). CINTEGERI 
(6) “XSHIFT’. THIS CARD ALLOWS THE USER THE OPTION OF GIVING 
THE X INPUT AS EITHER SHIFTED OR NORMAL. IF A SHIFTED X 
COORDINATE DESIRED THEN SPECIFY A ’1’ ON THIS CARD. IF 
NOT, GIVE A ‘0’. CINTEGERI 
(7) ‘YSHIFT’. THIS CARD ALLOWS THE USER THE OPTION OF INPUT- 
ING EITHER THE Y COORDINATE OR INCLINED SPAN DISTANCE 
ALONG THE SKEWED EDGE. IF SKEWED Y COORDINATE 2 DESIRED 
THEN SPECIFY A ’1’% ON THIS CARD. IF NOT, GIVE "O07. 
CINTEGER J 
(6) Ii CARD #5 IS ’0’% THEN GIVE ON #7 : 


OAAOAOGAOGAIAAGIOAIOAOAIOOHIaAAIAIANIANAANAANAONADA 


hook wk ako ok Oke oe oe OK OK OK oe oe oe ok ae ae 
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(Ad THE X COORDINATE OF THE WHEEL NUMBER GIVEN ON 
CARD # 4. CREALI 
IF CARD #5 IS ’1° THEN GIVE ON #7 : 
(B) THE X DISTANCE OF THE WHEEL FROM THE LEFT SKEWED 
EDCE OF THE DECK. CREALI 
(9) IF CARD #6 IS ’0’ THEN GIVE ON #8 : 
(A) THE Y COORDINATE OF THE WHEEL NUMBER GIVEN ON 
CARD # 4. CREALI 
IF CARD #6 IS ‘1° THEN GIVE ON #8 : 
(B) THE DISTANCE TO THE WHEEL AS MEASURED ALONG THE 
INCLINED 
SPAN EDGE OF THE DECK. CREALI 
HHENOTE### THE OPTIONS OF SKEWED COORDINATES ARE INDEPENDENT 
FOR THE X AND Y AXES. IT IS RECOMMENDED THAT THE SKEWED OPTION 
BE USED WHEN PLACING A WHEEL ON ANY GRID LINE OR MESH NODE, 
ALTHOUGH IT IS OFTEN CONVENIENT AT OTHER TIMES AS WELL. 
(10) ‘KEY’. IF KEY=0 THE COORDINATE MESH WILL NOT BE PRINTED 
OUT. IF KEY =1 IT WILL. CINTEGER] 
(11) “’NXSTEP’. THE NUMBER OF ELEMENTS IN THE X DIRECTION. 
CINTEGER I 
€i2) “NYSTEP’. THE NUMBER OF ELEMETS ALONG THE (SKEWED) Y 
DIRECTION. CINTEGER]J 
C13) °XINC’. THE X INCREMENT (INCHES). CREALI 
(14) ‘“YSKINC’. THE INCLINED Y INCREMENT (INCHES). CREALI 


o— oe see coe oe ees oo oo ee oe ee cee ee ee ce wee ec) ee oom SS et et er ce Ss ee ee Oe oem ee comme ee ec ee oe ee Se SE ee ee Se ee ee 


COICO OI Oe CCC CIC CCC Ie ie CIC et 


eS oes ees Oe oe oem See oem GP ee re ee cee ee ee ee ee ee ee ee ee ee ee ee ee CS a ee ee ee eee ee 


PLACE WHEEL # 3 ON NODE (5,4) OF A 30 DEG SKEWED DECK. THE 
DECK MEASURES 540" kK 340" (X DIMENSION AND SKEWED SPAN, RESP) 
AND THE MESH SIZE IS 34" IN THE X DIRECTION AND 45" IN THE ¥ 
DIRECTION. THUS THERE ARE 10 INCREMENTS ALONG THE X AXIS AND 8 
INCREMENTS ALONG THE SKEWED Y AXIS (PLEASE SEE SKETCH). SAY 
THE TRUCK TYPE IS FDOT SU-4. WE WOULD INPUT AS FOLLOWS: 


CeCe ee eee ee eee ee ee ee ee ee ee ae 


CARD INPUT COMMENTS 
em ees ee ees eee gree wees Se res eee ee crm ee ms ree ee ee cn wees eS cme we cones ey ee es es wees AS es Oe RE OES SOD GU SD eS SH a Sm Se Sm He OS ES eS es eS me See OS OS ER Somes One Cum ts 
1 FDOT SU-4 EXAMPLE OUTPUT LISTING TITLe 
2 1 TRUCK TYPE IS FDOT SU-4 
3 30.0 SKEW ANGLE IS 3O DEGREES 
4 1 LOAD CASE NUMBER (USER ’S 
OPTION) 
2 3 INPUT WHEEL IS #3 
6 1 THE XSHIFT IS ON 
7 1 THE YSHIFT IS ON 
8 162.0 NODE (5,4) IS 162" FROM THE 
LEFT EDGE OF THE DECK 
9 180.90 NODE (5.4) IS 4 SPACES AT 


45" EACH CINCLINED DISTANCE) 
FROM THE BASE OF THE DECK 


10 0 KEY=0 MESH WILL NOT BE 
PRINTED OUT 

11 10 NUMBER OF X STEPS [S 10 

12 8 NUMBER OF Y STEPS [5 8 

13 24.0 X INCREMENT IS 94” 

14 49.0 Y (SKEWED) INCREMENT I5 45” 


eH He mH A AK Ee OK KE IK Fe KO OK OR He HR 


JELEIEIE ETE SES EEE EE ASEH HHH HHHHHKHHHHKKHHREHHHERHHHREEHKEKKRERRKERKERHEEE 
DIMENSION XCOR(30, 30), YCOR(30, 30); REACT(30, 30), 
$ WHEELX(18), WHEELY(18),WT(18), AXLSP(8), AXLWT(9) 
INTEGER ene SERRA HAUT, TRISEE 
CHARACTER ALPHA 
OPEN (UNIT=1,FILE=’C.PAULILOADER. DAT’, STATUS=“OLD*) 
OPEN gee. fel OUTEUT ALIS S16 | 
CALL TITLE1 (ALPHA 
CALL INITL(WUX, WJY, PHI, NXSTEP, NXCROS, XINC, NYSTEP: NYCROS, 
$ YSKINC, LCN, KEY, XSHIFT, YSHIFT, WNUM, AXLSP, NA, NAXLSP, NW, 
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BADIM RICE TARDE DOV 
EP, XINC, NXCROS, XCOR, NY nA , 
$ NYCROS, YCOR, KEY, ERR) aia Boece aa 
CALL LOCATE (CWJIX, WY, W1X, WLY, PHI, WNUM, ERR, OFFWT, XSHIFT, 
$ YSHIFT, USERX, USERY, AXLSP, NA, NSXLSP, NW. WIDTH, WT. NACROS, 
$ NYCROS, REACT, WHEELX, WHEELY, DDX, DDYSK) 
CALL TITLE2 
CALL SEARCH(PHI, WT, REACT, XCOR, YCOR, WHEELX, WHEELY, NXSTEP, 
$ WNYSTEP, XINC, YSKINC, NXCROS, NYCROS, NW) 
CALL RESULT (NXCROS, NYCROS, REACT, OFFWT, XSHIFT, YSHIFT, 
Se $ Be ee Vex, WY) 


END 
SUBROUTINE INITLCWUX, WIY, PHI. NXSTEP, NXCROS, XINC, NYSTEP, 
$ NYCROS, YSKINC, LCN, KEY, XSHIFT. YSHIFT, WNUM, AXLSP, NA, 
$ NAXLSP, NW, WIDTH, TRTYPE, WT, DDX, DDYSK) 
CHER ERHHR HEHE HKEKRHHHHRHRRRHHK EHH EREEEREKHEER 
C SUB INITL READS IN THE INPUT DATA AND SETS VARIOUS INITIAL % 
C VALUES BASED ON THIS INFORMATION. % 
CHEER EH EHH HHE RHE RRHHHAHHKSHERHREEHRREERRREEHE SHER EEECERE 
DIMENSION AXLSP(8),WT(13), AXLWTC9) 
INTEGER XSHIFT, YSHIFT, WNUM, TRTYPE, ERR 
WRITE (6,310) 
310 : ne) See ee ee - 


ee Cue EY Sune om ee we SOP Cue Garo Cue GE fee Core cee come SY ey ee ee ees ey ey ey cep eee ee ee ce 


READ (1,315) TRIYPE 
FORMAT (111) 
Be CIRTYPE.NE.1)? GO TO 325 
WRITE (4) 320) 
320 FORMAT (7, 7%, “TRUCK TYPE’, 18X, “FDOT SU-47%) 
CALL TRIYP1(CAXLSP, NA, NAXLSP, NW WIDTH, WT) 
GO TO 340 
329 Meet YPe.NNE.9) GO TO 360 
CALL CUSTM1 (NA, NAXLSP, Ni) 
CALL. CUSTM2(AXLSP, NA, NAXLSP, NW, WIDTH, WT, AXLWT ) 
WRITE (6,330) NA, WIDTH 
330 FORMAT (7, 7X, “TRUCK TYPE’, 19%, ‘CUSTOM’, 7/, 14%; 
‘NUMBER CF AXLES’, 8X,12,//,14X, “TRUCK AXLE WIDTH’, 
6X, Poe 2) 
DO 340 J=1, NAXLSP 
WRITE (6,335) I, AXLSPC(TI) 
335 FORMAT (/,14X, “AXLE SPACE # 7%, 11, 8X, F6.2) 
340 CONTINUE 
DO 350 I=1,NA 
WRITE (6,345) 1, AXLWTC(I) 
349 FORMAT (/,14X, “AXLE WEIGHT # ’, 11. 7X, F6.2) 
350 CONTINUE 
360 READ (1,362) PHI 
362 FORMAT (1F10.1) 
WRITE (6,344) PHI 
364 FORMAT ¢/,7X, “ANGLE OF SKEW (DEG) ’,11X,F4.1) 
READ (1,366) LCN 
366 FORMAT (112) 
WRITE (6,368) LCN 
368 FORMAT (/,7X, ’LQAD CASE NUMBER’, 14X, I2) 
READ (1,370) WNUM 
370 FORMAT (112) 
WRITE (4,372) WNUM 
372 FORMAT (/, 7X, ‘INPUT WHEEL NUMBER’, 12x, 12) 
READ (1,374) XSHIFT 
374 FORMAT (111) 
READ (1,376) YSHIFT 
376 FORMAT (1iI1) 
READ (1,378) WUX 
378 FORMAT (1F7. 2) 
READ (1,380) WJY 
380 FORMAT (1F7.2) 
IF (XSHIFT.EQ@.0) GO TQ 383 


i) 
ps 
Ui 


+f 4% 


oD 





382 


383 
384 


386 
388 


389 
390 
392 
393 
394 


396 
398 


405 


410 
412 
414 
420 
430 
432 


434 


930 


+h th th 


be 


FORMAT (7, 7%. XIN 
X, “X PUT (SHIFT 2p 

Cama dae, ON) °, 10X,F 7.2) 
WRITE (6,384) WJX 

FORMAT (€/, 7X, “°X INPUT (SHIFT OFF)’, 9X, F7. 2) 

IF (YSHIFT.EQ@.0) GO TO 389 
WRITE (6,388) WJY 
FORMAT mononeN ANPUT <SHIFT ON) %, 10X,F7. 2) 


WRITE (6,390) WJY 
FORMAT (7, 7X, “Y INPUT (SHIFT OFF)’, 9X, F7. 2) 
Rene <1, 393) KEY 
FORMAT (111) 

IF (KEY.EQG.0) GO TO 396 

WRITE (6, 394) 
FORMAT (7, 7X, “MESH PRINT REQUEST’, 13X, ’YES”) 
GO TO 400 
WRITE (6, 398) 
FORMAT (7/7, 7X, “MESH PRINT REQUEST’, 12X; ‘NO’) 
WRITE (6, 405) 


FORMAT S11, 7%) (errr er rr ee + 


READ (1,410) NXSTEP 
FORMAT (112) 
READ (1,412) NYSTEP 
FORMAT (112) 
WRITE (6,414) NXSTEP, NYSTEP 
FORMAT (7/7, 7X, “NUMBER OF ELEMENTS’, 9X, I2, 10X, 12) 
READ (1,420) XINC 
FORMAT (C1F6.2) 
READ ¢1,430) YSKINC 
FORMAT (1F6. 2) 
WRITE (6,432) XINC, YSKXINC 
FORMAT (7, 7X, “ELEMENT SIZE’, 14X,F6.2, 6X, Fé. 2) 
DDX=NASTEP*XINC 
DDYSK=NYSTEP*EYSKINC 
WRITE (6,434) DDX, DDYSK 
FORMAT (¢7, 7X, “DECK DIMENSION’, 11X,F 7.2, OX, F7. 2) 
NXCROS=NXSTEP +1 
NYCROS=NYSTEP+1 
RETURN 
END 


= 


= 


95 


SUBROUTINE GMESH(PHI, NXSTEP, AINC, NXCROS, XCOR, NYSTEP, 


YSAINC, NYCROS, YCOR, KEY, ERR) 


CHHERHHHHHEHHRRHHHHKHHRHEPHHKHHKHKRHHHRHHKRHRHEEHEHHRHKRRERRERERKRRKEREHREK IR 


C SUB GMESH GENERATES THE SKEWED MESH AND ALLOWS FOR 


C PRINTOUT FOR KEY=1. 
CHEKHHHAKHAHHHE KE HRHHKHEKRHR HHH HHKHHHRHHHERRKHEHE EE RHEER EEE RHRHREREEE HE 


750 
270 


DIMENSION XCOR(NYCROS, NXCROS), YCOR(NYCROS, NXCROS) 


INTEGER ERR 

PHI=3. 14159265359#PHI/1i80. 
X=0. 0 

Y=0.0 

A=0.0 

XSTEP=0.0 

NCOUNT=0 

ERR=0 

DO 9970 J=1,NXCROS 

DO S80 I=1,NYCROS 
NCOUNT=NCOUNT+1 

YCOR (I, Jd=¥Y 

XCOR CI, J) =X 

IF (KEY.EQ@.0) GO TO 570 
WRITE(6, 560) NCOUNT, I, J, X,Y 


FORMAT (/,5X, ‘NODE # %,13,5X, CI, J)=€%, 12, %, %. 12, 


ye pie 7 wes OX> “Y=",F7.2 
A=A+YSKINC 


7) 


OP TIGNAL 


“y DX» 


* 





980 


970 


940 


96 


X1=A#SINC(P 
¥1=A*COS(P 
X=XSTEP~X1 
Y=Y1 

CONTINUE 
A=O., O 


¥=0.0 

XSTEP=XSTEP+XINC 

X=XSTEP 
CONTINUE 

RETURN 

END 
SUBROUTINE LOCATE (WUX, WY, WIX, W1LY, PHI. WNUM, ERR, OFFWT, 


$ XKSHIFT, YSHIFT, USERX., USERY, AXLSP, NA, NAXLGP, NW. WIDTH, WT, 
$ NACROS, NYCROS, REACT. WHEELA, WHEELY, DDX; DDYSK) 


HI) 
HI) 


CHHHMHHHHHHHHEHHHHHHHHHHE ERR ERHHHHHHHREEEHE HICH PE RHEEHKKHERHHREE EEE 
SUB LOCATE USES THE INPUT WHEEL COORDINATES OF WJUX AND WJY TO * 
FIND THE COORDINATES OF THE LEFT FRONT (i.e... DRIVER’S) WHEEL. * 
THEN IT WILL LOCATE THE CQGORDINATES GF THE REST BASED ON THIS # 
POSITION. IT ALSO CHECKS TO BE CERTAIN THAT EACH WHEEL FALLS # 


os 
TRE DECK ITSELF AND PRINTS A WARNING LIST FOR THOSE WHEELS # 
¥ 


THAT ARE NOT WITHIN THE DECK BOUNDARIES. THE TOTAL WEICHT OF : 


THE WHEELS THAT ARE OFF THE DECK WILL BE CONSIDERED IN THE + 
~~ 


STATICS CHECK AT THE END OF THE PROGRAM. 
HHEHKHKKEKHHKKAHHKHHHHKRRHHEHHEHEERRHRHHKHHRHEEEEE RH EHRHEHKRHHRHHEEREEE 


C 
C 
C 
c 
C CN 
C 
. 
c 
C 


600 


602 
507 
910 
942 


944 


946 


948 


990 


992 


DIMENSION WTCNW). AXLSP(NAXLSP), WHEELX (NW), WHEELY (NW), 
REACT (CNYCROS, NXCROS) 

INTEGER WNUM, ERR. XSHIFT, YSHIFT 

IF (WNUM.LE.NW) GO TO 602 

WRITE (1,600) 
FORMAT (//,. 7X, ’### ERROR ##*® ILLEGAL WHEEL NUMBER INPUT’ > 
ERR=1 

GO TO 1070 
USER X=WUK 

USER Y=WJY 

ae Cronae Te EQ. 0) CO TO 607 

WIY=WJIYXCOS (PHI ) 

IF (XSHIFT.EG.0O) GO TO 710 
WUX=WIX-(WIY#SINCPHI) /SINC1. 570796327-PHI )) 
Z1=WIDTH#COS(PHI ) 
Z2=WIDTHS#SINC(PHI) 

IF (WNUM.NE.1) GO TO 942 
WI XK=WUX 

WIY=WUY 

GO TO 980 

IF (WNUM.NE.2) GO TO 944 
W1A=WUX—-Z1 

WLY=WJUY-Z2 

GO TO 980 
A=AXLSP (1) 

IF (WNUM.NE.3) GO TO 946 
WIX=WUX-CA)#SINCOPHI) 
WIY=WJUY+ (CA) #COS(PHI) 

GO TO 980 

IF (WNUM.NE.4) GO TO 948 
WIX=WIUX-C CAD BSINCOPHI)D+Z1) 
WIiY=WJUY+( (A) #COS(PHI)-Z2) 
GO TO 980 
B=AXLSP (2) 

IF (WNUM.NE.5) GO TO 950 
W1X=WJIX-CAt+B)#SINCPHT) 
W1Y=WJY+ (A+B) #COS(PHI) 
GO TO 93s0 


™~ 


IF (WNUM.NE.&5) GO TO 952 
WiX=WIX-(CA+BI#SINCPHI)+Z1) 
WIY=WJY+( (A+B) #COS(PHI)-Z2) 


GO TO 780 
C=AXLSF (3) 
IF (WNUM.NE.7) GO TO 954 





954 


9956 


998 


950 


762 


764 


966 


768 


970 


972 


974 


976 
978 


730 


oy 


Wi X=WJUX-(A+B+C ) #SIN( PHT 

LY=WJUY+(A+B+C ) #COS (PHI 

O TO 980 
id (WNUM.NE.8) GO TO 9565 
1TX=WUX-¢€ (At+B+C ) #SIN(PHI)+Z1) 
LY=WJUY+¢(At+Bt+C) #COS (PHI) -Z2) 
GO TQ 380 
=AXLSP (4) 
-F (WNUM. NE. 9) GO TO 3758 
1X=WUX-(CAt+BtC+D) #SINCPHI) 
1Y=WJY+ (A+B+C+D)#COS(PHI) 
GO Te 380 

IF (WNUM. NE. 10) GO TO 940 
Wi X=WJUXK-( €44+B4+C+D) #5 IN¢(PHI ) 
W1IY=WJY+( €44+B4+C+D) #COS (PHI) 
GO TO 780 

E=A4XLSPC(5) 

IF (WNUM.NE.11) GO TO 3762 
WiX=WUX-(CAt+Bt+C+D+E) X¥SIN(P 
WLIY=WUY+ CAt+Bt+C+Dt+E) #COS(P 
CGO TO 980 

IF (WNUM.NE.12) GO TO 9764 
Wix=WUX—-¢( C(At+B4+C4+D+E) #SINCPHI 
WIY=WJIY+¢ CA+Bt+C4+D+E) #COS (PHI 
GO TO 980 

F=AXLSP (6) 

IF (WNUM. NE. 13) GO TO 766 
WIX=WJUX-CAt+Bt+C+D+Et+F ) #SINCPH 
WLY=WJUY+( A+B+C+D+E+F ) #COS(PH 
GO TQ 980 

IF (WNUM. NE. 14) GO TO 968 
WIixX=WJUX-C CAtB+C+D+E+F ) #5I 
WIY=WJIY+( (At+Bt+C+D+E+F )#CO 
GO TQ 980 

G=AXLSP (7) 

IF (WNUM.NE.15) GO TO 970 
WiX=WIUX-(At+B+C+D+E+F+C) #SIN(P 
WILY=WJY+ (4+B4+C+Dt+Et+F+¢6)#COS(P 
CO TA ¥Y8BO 

IF (WNUM.NE.16) GO TG 972 
WiX=WUX-( CA+B+C+Dt+E+F+6) #SIN( PHI) +213 
WIY=WJUY+( (4+B+C+D+E+F+G6)#COS(PHI)-Z2) 
GO TO Lae 

H=AXLSP (8) 

ie CHINUMS NE.17) GO TO 974 
WixXx=WJUX-(CA+Bt+C+D+E+F+6G+H) #SIN(PHI) 
W1Y=WJUY+ (A+B+C+Dt+E+F+6+H) XCOSC(PHI ) 
GO TO 980 

IF (WNUM. NE. 18) GO TC 976 

Wi X=WJUX-( (At+B+C+D+E+F+G6+H) SINC PHI 
WiY= pee COTE SES Ot) COSKE HS 
GO TO 980 
WRITE (20978) ,, 
FORMAT OI a All ILLEGAL WHEEL NUMBER INPUT 
‘ ON DATA CARD # 107 


ERR=1 
GO TQ 1090 


W 
GC 
= 
W 
W 
GC 
D 
I 
W 
W 


wet ) 
~Z2) 


gee 


ez 


d+Z 
a7. 


1) 
=o) 


=O 

WHEELXC1)I=W1X 
WHEEL XK (2)=W1X+Z1 
WHEELYCLII=HWLY 
WHEELY (2) =W1Y+Z2 

T={[+2 

ir Wi. EQ. NW) GO TO 1010 
A=AXLSP (1) 

WHEELX (3) =W1X+(A)#SIN(PHI) 
WHEEL X (4) =WHEELX(3)+Z1 

WHEELY (3) =W1Y-AXCOSCPHI) 
WHEELY (4) =WHEELY (3)+Z2 





990 


98 


T=I+2 

IF (I.EQ.NW) GO TO 1010 
LG=AXLSP (2) 

WHEEL X (CD) =W1X+( A+B) #SIN(PHI ) 
WHEEL X(6)=WHEELX(5)+Z14 

WHEEL Y (S)=WiY-(A+B)+COS(PHI ) 
pee eC EELY (5) +22 

oo 

IF (I.EQ.NW) GO TO 1010 
C=AXLSP (3) 

WHEEL X (7) =W1IX+(A+B+C)#S INC PHI ) 


WHEEL X (8) =WHEELX(7)+Z1 

WHEELY (7)=W1Y-(4+4+8+C )#COS( PHI) 
pee CY 7) +22 

= + 

IF (I.EQ.NW) GO TO 1010 

D=AXLSP (4) 

WHEEL X(9)=W1kK+(4+B+C+D)#SIN(PHI ) 
WHEEL X C10) =WHEELX(9)+Z1 

WHEELY (9)=W1Y-(A+B+C+D)#COS (PHI ) 
eee oe) eee 

=JT+ 

IF (€I.-EQ@.NW) GO TO i010 

EF=AXLSP (5) 

WHEEL X(11)=W1X+(A4+B4+C+D+E) #SIN(PHI ) 
WHEELX (12)=WHEELX(11)+2Z1 
WHEELY (11)=W1Y—-(At+B4+C4+D+E) X#COSCPHI ) 
WHEELY( 12)=WHEELY(11)+Z2 

l=J[+2 

IF (I.EQ.NW) CGO TO 1010 
rF=AXLSP (64) 

HEEL XC LIS)I=Wik+CAtE+C+D+E+F DXSINCOPHT ) 
WHEEL X (14) =WHEELX(13)+Zi 

WHEELY (13) =W1Y-(ATBTtC+D+E+F J X#COS(PHI) 
as ae YS S22 

= [+ 

IF (IJI.EQ.NW) GO TO 1010 
G=AXLSP (7) 
WHEEL X (15)=W1X+(A+B+C+D+E+F+G) #SINOPHE D 
WHEEL X (C16) =WHEELX(15)+Z1 
WHEELY(15)=W1Y- (A+B+C+D+E+E+C) #COS (PHI D 
eae YS to) re 

=[+ 

IF (I.EQ.NW) GO TO 1010 

H=AXLSP (8) 
WHEEL X (17) =WiX+(At+B+C4+Dt+E+F+CGt+H) #SINCPHE ) 
WHEELX (16) =WHEELX(17)+2Z1 
WHEELY(17)=W1iY- (A+B+C+D+E+F+6+H) #COS(PHE ) 
eee ICEL E7) +22 

=J+ 

IF (I.EQ@.NW) GO TO 1010 

WRITE (6,990) 


FORMAT (//,7X, ‘ILLEGAL NUMBER OF WHEELS INPUT’) 
GO TO 1090 


1010 eS ae =0.0 


1080 I=1,8 
A=WHEELY(T) 
IF (A.GE.0.0) GO TO 1015 
OFFWT=OFFWT+WT (CT) 
WTCI)=0.0 
RITE (46,1013) I 


W 
5 ae FORMAT (///, 7%, (HR#HEWARNINGH SHH WHEEL # %, 12, 
I 


S OFF THE DECK TO THE SOUTH’) 
GO TO 


10380 
1015 Se ae 


(A.-LE.-PSD) GO TO 1025 
OFFWT=OFFWT+WT (CT) 
WT(1)=0.0 
WRITE (46,1020) I 





1020 
$ 


10295 


19330 
$ 


1035 


1040 
$ 


1080 
640 


690 
1090 


$ 


gag 


FORMAT (///, 7X, “####WARNING#HR## WHEEL # 7’, 12, 
“ IS OFF THE DECK TO THE NORTH’) 
GO TO 1080 
OF FSET=AXSIN(PHI)/SIN(1.570796327-—PH! ) 
B=WHEELX (I) 
B1=DDX-OFFSET 
IF (B.LE.B1) GO TO 1035 
OF FWT=OFFWT+WT (I) 
WrCI)=0.0 
WRITE (6,1030) I 
FORMAT (///, 7X, (##R#WARNINGH### WHEEL # °, 12, 
“IS OFF THE DECK TO THE EAST’) 
GO TQ 1080 
B1=0. O-OFFSET 
IF (B.GE.B1) GO TO 1080 
OFFWT=OFFWT+WT (TI) 
WrCI)=0.0 
WRITE (6, 1040) I 
FORMAT (///, 7X. (####WARNINGHEHX WHEEL # %, I2, 
eeloeerts THE DECK TO THE WEST‘) 
CONTINUE 
DO 690 I=1, NYCROS 
DO 650 J=1,NXCROS 
REACT(I, J)=0.0 
CONT INUVE 
RETURN 
END 
OUBROUTINE SEARCH(PHI, WT, REACT, XCGR, YCOR, WHEELX, 
WHEELY, NASTEP, NYSTEP, XINC, YSKINC, NACROS, NYCROS, NW) 


CHEHHHHHAHHKHKRHHEHHHEKRHKEHRKEKKHHKHKHHKKRHHHKKHSHHKKRHHEHHREHKHKHKEHHRREREE 


meow S 
C RESPE 
© VERTI 
C OF TH 
C SEARC 
C DISTA 
C NODES 
= mer D 


EGRCH IS USED TO FIND THE LOCATION OF EACH WHEEL WITH 
CT TO ITS SURROUNDING NODES. DO THIS BY FIRST SEARCHING 


CALLY TO GET DISTANCES Di & D2 WHICH ARE THE COQRDINATES 


FEF MODES IMMEDIATELY ABOVE & BELOW THE WHEEL, RESP. THEN 
H HORIZONTALLY ALONG THE Y COORDINATE OF D2 TO GET THE 


NCES D3 AND D4 WHICH ARE THE PROPER XK COORDINATES OF THE 


ALONG THIS LOWER ’TRACK’. THEN USE A SIMILAR PROCEDURE 
Y AND D& ALONG THE UPPER ‘TRACK’. 


KK OK HK MK OK OK 


HHERHHHHKEEHHKHHHHHHHRHHERKRRHERHHHHHEHHPHERHRHERHRRRHKHEEREHKRREPHEHKRHKERREEER 


$ 
$ 


705 


DIMENSION WT (NW), REACT (NYCROS, NXCROS), 
XCOR (NYGROS, NXCROS), YCOR(NYCROS, NXCROS), WHEEL X (NW), 
WHEELY (NW) 
DO 850 K=1, NW 
WICHK=WT (K) 
IF (WTCHK.EG.0.0) GO TO 850 
J 
LCOUNT=0 
C=YCOR(T, 1) 
D=WHEELY (K) 
CALL RNDOFF(C, D) 
IF (C.EQ@.D) GO TO 710 
IF (€.6T.D) GO TO 720 
+] 
LCOUNT=L.COUNT +4 
co To 
aT 
i 
D1=YCOR(I1, 1) 
eo0To) 721 
it 
I2=11-1 
NNE=T1 
NSE={2 
Da=YCOR (11, 1) 
D2=YCOR(I2, 1) 
yoo TO 729 
LCOUNT=LCOUNT-NYSTEP 
F=WHEEL X (K) 





722 


723 


724 


729 


725 


729 
730 


731 


732 


$ 


100 


G=XCOR(T1,N) 
CHECK3=ABS (WHEEL X (K)-—(N-1)#XIN 
CHECK4S=ABS(WHEELY (K)#SIN(CPHI)/ 
CALL RNDOFF (CHECKS, CHECK4) 
CALL RNDOFF(G, F) 

IF (CHECKS. EQ. CHECK4) GO TO 723 
IF (G.CGT.F) GO TO 725 

N=N-+ 1 

L_COUNT=LCOUNT+NYSTEP 
NE=NE+NYCROS 

GO TO 722 

TS=N 

DS=XCORC(I1, 15) 

W=WT CK) 

A=WHEELX (CK) 

B=WHEELY (K > 

WRITE (6,724) K,W, LCOUNT, NE, A,B 

FORMAT (/7,9X, 11, 6X,F5.2., 4X, ‘ABOVE ’,12,8X, “OQNLY /%, 

IZ, 7X, F7.2, 2X, F7.2) 

RS=WT (CK) 

REACT(CII, IS)=REACTCII, 15)+RS 
€0 TO 850 


C) 
SINC?. 970796327 -FH1)) 


0 
A=WHEELX (CK) 
B=WHEELY (K) 
WRITE (6,726) K, W, LCOUNT, NW, NE, A, B 
FORMAT (7,9X, 11, 6X,F5.2, 4X, “ABOVE %,I2,8X, “ONLY %, 
I2, 1X, 12, 4X,F7.2) 2X, F7.2) 

RL=(DS-A) FW7 (CK) / CDS—-DSG) 

RR=WTCK)-RL 

REACTCI1, IS)=REACTCI1I, I5)+RR 

REACT(I1, I6)=REACT( TI, 16)+RL 


a TO 850 
LCOUNT=LCOUNT-NYSTEP > 


E=XCOR(T2, L) 

F=WHEELX (CK) 

CHECK4=ABS (WHEELY (KO #SINC(PH 
CHECKS=ABS (WHEELX (CK) -(CL—-1)% 
CALL RNDOFF (CHECKS, CHECK4) 
CALL RNDOFF (CE, F > 

IF (CHECKS.EQ. CHECK4) CGO TO 73i1 
IF (E.GT.F) GO TO 740 

L=L+1 

LCOUNT=LCOUNT+NYSTEP 
NSE=NSE+NYCROS 

GO TO 730 
I3S=L 

IS=13 
D3=XCOR(I2, J 
DS=XCORC(I1, 1 
NB=NSE 
NA=NB+1 
W=WT CK) 
A=WHEELX (CK 
B=WHEELY (K 
WRITE (4, 7 
FORMAT (¢/, 


I>) /SINC(1. S7O0796327-PHI) ) 
XINC) 


3) 
3) 


F5_2,4X, ‘RT OF ’,12,8X, “ONLY 7’, 
- 


I 





740 


790 
750 


770 


780 


790 
F929 


797 
$ 


10] 


XOFF1=(D-D2)+#S 
CHECK1=XCOR(I1 
WXOFFLI=CHECKi 
CALL RNDOFF (W 
IF (WXOFF1. GT. 
ISG=L+1 
I4=L 
LCOUNT=LCOUNT+NYSTEP 
NSE=NSE+NYCROS 
GO TO 7&0 

I3=L 


T4=L-1 
D3=XCOR(I2, 13) 
D4=KCOR(I2, 14) 
NSW=NSE-NYCROS 


N=1 
G=XCOR(I1,N) 

CALL RNDOFF(G, F) 

IF (G.GT.F)} GO TO 780 

N=N+1 

NNE=NNE+NYCROS 

GO TO 770 
XOFF2=(D1-D)*SIN( PHI) /SIN(1.570796327-PHI1) 
CHECK2=XCOR(T1,N) 

WXOFF2=F—-(CHECK2—-XINC ) 

CALL RNDOFF(WXOFF2, XOFF2) 

IF (WXOFF2.CT.XQFF2) GO TO 790 

I9=N-1 

I6=N-2 

NNE=NNE-NYCROS 

LCOUNT=LCOUNT-NYSTEP 

GO TO 795 
TO=N 
I6=N-1 
DO=XCORCI1, 
D6=KXCORCI1 
NNW=NNE-NYC 
W=WT (CK) 
A=WHEELX (CK ) 

B=WHEELY (K ) 

WRITE (6,797) K, W, LCOUNT, NSW, NNW, NSE. NNE, A, B 
FORMAT (¢/7,9X, 11, 6X, F5.2, 7X, 12,10X,12,1X, 12, 1X, 12, 
1X, 12, 4X,F7.2,2X,F7.2) 


(SHI) /SIN(1.570796327-PHI) 
FF 
OF 


t 
2) 
x 1, 
ae 


IN 
- 
0 XOFF 1) 

A } GO TO 750 


2) 
t6 
RO 


) 
= 


CHEK HHHHERERHEERHEREH AH HAAHAKHHHHEEK REE HERR RKKEKEREKHHERKKEEH RRR 
COMPUTE NODAL POINT REACTIONS BY PUTTING AN IMAGINARY STRINGER * 


C 
C 
C 
C 
C 
C 
C 
C 


HORIZ 
STATI 
OTRIN 


LOADS. 
NG TOTAL OF EACH WHEEL ’S CONTRIBUTION IN A STORAGE MATRIX 


RUNNI 
REACT 


HHEEKRHRHHHHHHKHKKHKHRKKKKHRKKKKRKHKHHHKRKHHHKRHKEHHHKERAHRHHHKKHHERRERHEEHERE 


BOO 


890 


9135 


ONTALLY ACROSS THE ELEMENT AND FINDING THE EQUIVALENT 
C LOADS ON THE SKEWED EDGES. THEN PUT TWO IMAGINARY 


GERS ALONG THE SKEWED EDGES AND SOLVE FOR THE NODAL Bs! 


AS THE REACTIONS FROM EACH WHEEL ARE FOUND, KEEP 
(I,J). 


D7=WHEELX (K)—-D4 

DS8=(WHEELY (K)-D2) #SIN(PHI)D/SIN(1. 570796327-PHI1 ) 
RL=WT(K)#¢ (D3—-D4)—-(D7+D8) )/ (D3-D4) 

RR=WT CKI-RL 


RNNW=RL# CWHEELY (K)-D2)/(D1—-D2) 
RNSW=RL#(DI-WHEELY(K))/(D1-D2) 
RNNE=RR# (CWHEELY(K)-D2)/(D1i-D2) 
RNSE=RR*(DI-WHEELY(K))/¢(D1-D2) 
REACT(I1, I6)=REACTCI1, 16)+RNNY 
REACT (I2, 14)=REACT (12, 14)+RNSW 
REACT(I1, I5)=REACTCI1, 15) +RNNE 
REACT(I2, 13)=REACT(I2, 13) +RNSE 
CONT INVE 

RETURN 


END 
SUBROUTINE RNDOFF (X,Y) 


* 


* 
= 
* 
os 
% 
*% 


CEH HKRHHHKKKRHHESREHHEHH RHE RHHKHKHHHHHHRKKHKHRHEHHERHKRERHKHKRRHEREKRREERS 
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C SUBROUTINE RNDOFF ROUNDS OFF THE VALUES ASSIGNED TO X AND Y t% 
© PO THREE DECIMAL PLACES. + 
HEHRKKERHRHHRHHRHHHHHEER HERR HKRHHEHHHEHEHRH EEE HEKHRHERHHKHREEHHERHEK 

A=1000. 0#X 

Y=1000. O#Y 

X=ANINT(X)71000. 0 

YH=ANINT CY) /1006. 0 

RETURN 

END 

SUBROUTINE RESULT (NXCROS, NYCROS,; REACT, OFFWT,. XSHIFT- 

$ YSHIFT, WNUM, USERX, USERY, WUX, WY > 

CHEHHRHHKSHKHHHHHHHHHRESHEHHHHHRHHHHKEHHHHHEHEHHRHKHHERHHKERERHE HEHEHE 
C SUB RESULT OUTPUTS THE RESULTANT (NON-ZERG) NODAL LOADS AND % 
C CHECKS BY STATICS TO ASSURE THAT THE SUM OF THE WEIGHTS OF % 
C THE ’OFF-WHEELS’ ADDS UP TO THE TRUCK WEIGHT. * 
CHEHHHRHHHHAHHHHHHHHHE KEK RHR HHKHHHHERRHRHHHHHEHKERRHRHKHH RRR HEREERHERE 

DIMENSION REACT (NYCROS, NXCROS) 

INTEGER XSHIFT, YSHIFT, WNUM 

TOTAL=0. 0 

NCOUNT=1 

WRITE (6,890) 
B29 FORMAT (717%, ////;, 7%, “RESULTANT NODAL LOADS’, //,7X; “NODE #7’, 

$ OX, “LOAD (KIPS)’%,/) 
DO 880 J=1,NXCROS 
DO 870 I=1,NYCROS 


WRITE (6,897) NCOUNT,R 
897 FORMAT (/, 8X, 12, 9%, F7.2) 
360 NCOUNT=NCOUNT +1 
370 CONTINUE 
330 CONT INVE 
TRAWT=TOTAL+OFFWT 
WRITE (6,890) TOTAL, TRAWT 
BIO FORMAT (////, 7X, “STATICS CHECK’, //, 7%; 
‘TOTAL LOAD ON DECK = %,F9.2, °% KIPS’s 4/4; 7X, 
‘TRUCK WEIGHT = ’%,F5.2, ’ KIPS’) 
IF (XSHIFT.EQ@.0O) GO TO 900 
WRITE (6,899) USERX, WJUX 
BIS FORMAT (//, 7X, (####NOTERH##’, //, 7X, “SINCE THE XSHIFT’, 
‘ IS ON X INPUT OF’,F7.2, ’ GIVES X COORD /’,F7.2) 
200 feo yonirt-£G.0) CO TO 910 
WRITE (6,905) USERY, WY 
909 FORMAT (//, 7%, (‘H#HHNOTEX###’, //; 7X, “SINCE THE _YSHIFT’, 
$6 ’ IS ON Y INPUT OF’,F7.2,’ GIVES Y COORD ’,F7.2) 
710 SE ee 


E 

SUBROUTINE TRTYP1¢(AXLSP, NA, NAXLSP, NW, WIDTH, WT) 
CHHEHEHHH EEK KHEHHHHHHEHEHEHHHHHHHHHHHHHHHEHHHHHHHHREERAHHERHRERARHER REE ER 
C SUB TRTYP1 HOLDS THE REQUIRED INPUT DATA FOR THE FDOT SU-4 % 
C TYPE TRUCK cS 
CHEHHHHHHEKHEHEHHHHE HEHEHE KHER KEHHHHHEHHRERHERHE RHE RRR EER ERE 

DIMENSION AXLSP(3),WT(8) , 


NA=4 

NAXLSP=3 

NW=8 
WIDTH=76. 0 
AXLSP(1)=110.0 
AXLSP (2)=90.0 
AXLSP(3)=50. 0 


Th 


ti} 


9090 CONTINVE 
RETURN 


END 
SUBROUTINE CUSTM1 (CNA, NAXLSP, NY) 





HS 


CHER EERIE EEE EEE HIER EEE IEEE IESE SEE EERE ERK EES EREE 
C SUB CUSTM1 READS IN THE NUMBER OF AXLES AND FINDS THE NUMBER # 
C OF AXLE SPACES AND THE NUMBER OF WHEELS. THE INFORMATION IS + 
C REQUIRED FOR THE VARIABLE DIMENSIONS IN SUB CUSTM2. % 
CHEREHHAR RHR RRR IE ARIE ESCH LEASE RRHEEHEEHEE ESE ESCHER 

READ (1,705) NA 
705 FORMAT (111) 

NAXLSP=NA- 1 

NW=2+NA 

RETURN 

END 


SUBROUTINE CUSTM2(AXLSP, NA, NAXLSP, NW, WIDTH, WT, AXLWT) 
a ge, RARE SARAH RHE SRR a Rie eae ee 
REST OF THE REQUIRE zi ! + 
C TYPE OF TRUCK ae eames 
CRRA EERE IE SESE EEE KEE IEEE EH EEER IE EHLIEHIEHERERREKEKE 
DIMENSION AXLSP (NAXLSP), WTCNW), AXLWTCNSD) 
READ (1,740) WIDTH 
740 FORMAT C1FS. 2) 
DO 770 I=1, NAXLSP 
READ ¢€1,750) AXLSPCTI) 
790 FORMAT (C1F6. 2) 
770 CONT INVE 
DO 780 I=1,NA 
READ (1,775) AXLWT(C1) 
oI FORMAT (C1F6. 2) 
780 CONTINUE 
DO 790 K=1,NA 
I=2#K—-]j 
JEH2ZEK 
WT CT I=AXLWT CK) /2 
WTC JI=WECT) 
790 CONTINVE 
RETURN 
END 
SUBROUTINE TITLE] CALPHA) 
CHE HHHKHHERKHEHHRERREEHE SHE KHEHHHHHKRKRHE HH EHR HER HEHHRERHHKRKESHKSHRESEHEE 
. cra READS IN AND ECHO PRINTS THE TITLE GF THE OUTPUT + 
+ 
CHK HHI RHHHRKKAHERKAKHHHHRHHEHHKHHHEKHRERHEKKRHKHHAR HARE HREHERE 
CHARACTER ALPHA*77 - 
READ ¢1,200) ALPHA 
200 FORMAT (¢A77) 
WRITE (6,210) ALPHA 
210 FORMAT ¢€//7/7,7%X,A7/7) 
RETURN 
END 
SUBROUTINE TITLE2 
WRITE ¢€6, 700) 
700 FORMAT ©4775 DOXs SHREK RHE HHRHHHKHEHEHHEHHHER O, /, POX, 
FHM SESE KEE IEEE EEE HHHHHHHH SY, SJ, DOX, “#7, ABX. °H# 7%, /, 20X;, 
‘% RESULTS OF LATEST ANALYSIS #7. /, 20X, *’#®’%,28X, 7%’, /, 2OX> 
PRI EI RHHRHHHREHHHHEHH LY, Sf, POX, 
SHIEH Ee ee eee eH, S//, 7X, “WHEEL #7, 
2X, ‘WEIGHT’, 2X, “IN ELEMENT #%, 2X, “AFFECTING NODES’, 2X, 
“X COORD’, 2X, “¥Y COORD’, /) 
RETURN 
END 


aft tf 7 4 UF oF 
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CHHEHHRHEHHHERHKRKRKEHKKHEAKEHEERARKHKHHHHKHKRHRKRAEREEHEEHHKHKRKHREHRHKHE HHH HRRE 
C 


CREATE. FOR +t 


CHEB HAHHHRHEHKHRRREHHHRAMHHCRKRARRHRHHKHHAHHKKRHK FHHHHKHRRHRHRHRHEHHEHE EERE 


30 


tt ut 


CHARACTER BUFFER#*9, BETA#79, ALPHA#=44, GAMMA#15 

OPEN (UNIT=1, FILE=’CL.PAULIOQUTPUT. LIS’, STATUS=’OLD’) 
OPEN (UNIT=2, FILE=’C.PAULIPS20100. DAT’, STATUS=’OLD‘%) 
OPEN (UNIT=3, FILE=’C.PAULIPS2ZO100. DAT’, STATUS= "NEW ) 
DO 10 I=1, 1000 

READ (1,5) BUFFER 

FORMAT (3X; A?) 

IF (BUFFER.EQ. “RESULTANT’) GO TQ 15 

CONTINUE 

READ (1,20) Nw, FU 

FORMAT (////, 3%, 13, 11X%,F9.2) 

FU=—-FJ 

READ (2, ’€A79)’) BETA 

WRITE (3, °(A79)%) BETA 

READ (2, ’(A44, 1F5.1,A15)°) ALPHA, POSIT, GAMMA 
POSIT=POSIT+t1i2.0 

WRITE (3, ’(A44,F5.1,A15)°) ALPHA, POSIT, GAMMA 

DO 30 I=1,23 

READ (2, ’(A79)°) BETA 

WRITE (3, °(A79)°) BETA 

CONTINUE 

WRITE (3,40) NJ, FJ 

DO SO I=1,25 

READ (1,35) NJ FJ 

FORMAT (7/7, 3X, 13, 11X,F5. 2) 

FU=—-F JU 

IF--{NJ.EQ.0) GO TO 70 

WRITE (3,40) NJ, FU 

FORMAT (13, 1X. “FORCE Z ’,F0.2) 

CONTINUE 

WRITE (3,80) 

FORMAT (’STIFFNESS ANALYSIS’, /, 

‘CALCULATE AVERAGE PRINCIPAL BENDING RESULTANTS’, /, 
‘FINISH %} 

STOP 

END 





105 


ee 2 ee ee ce ee ee ce ee ee ee ee cee ere cme ae ree ce cee ee eee cee cme Cem cee cee ee ee Ow 


= a) Coe ee ee oo cee Oe Ce oe oe ome ce oe eee com cme cee wre ce cee oe ee: oe — 
-— — — eS a CS SUD CUED Se Oe cme cee cee cee ee cee cme ere cee SD cme SD “ee Sree ee cme cme GD ees Cee ee ee ee cee oe De eee ee ee 


$ GISTRUDL C.PAULIPS20100. DAT; PR=C£.PAULIPS201001.LIS/N 
$ @C.PAULIRUN 

$ GTSTRUDL_ C.PAULIPS20100. DAT; PR=LE.PAULIPS201002.LIS/N 
$ @L.PAULIRUN 

S GTSTRUDL £.PAULIPS2O100. DAT; PR=L£.PAULIPS2ZOLOO3.LIS/N 
6 @L.PAULIRUN 

$ GTSTRUDL C.PAULIPS20100.DAT; PR=C.PAULIPS201004.LIS/N 
S @L.PAULIRUN 

$ GCTSTRUDL [C.PAULIPS20100. DAT; PR=CL.PAULIPS201005.LIS/N 
$ @L.PAULIRUN 

$ GTSTRUDL C.PAVULIPS2Z0100. DAT; PR=L.PAULIPS201006. LIS/N 
$ @L.PAULIRUN 

= GTSTRUDL C. PAULIPS20100. DATi PR=C.PAULIPS201007.LIS/N 
$ GTSTRUDL C.PAVULIPS20100. DATi PR=E£L.PAULIPS201008. LIS/N 
$ @£.PAULIRUN 

$ GTSTRUDL C.PAULIPS20100. DAT; PR=L.PAULIPS201009.LIS/N 
§ @€C£.PAULIRUN 

$6 GTSTRUDL [.PAULIPS20100. DAT: PR=L.PAULIPS2010010.LIS/N 
$ @L.PAULIRUN 

> GTSTRUDL C.PAULIPS20100. DAT; PR=L.PAULIPS2010011.LIS/N 
$ @L.PAULIRUN 

$ GTSTRUDL [.PAVULIPS20100.DATi PR=£.PAULIPS2010012.LIS/N 
S$ @C.PAULIRUN 

$ GTSTRUDL C.PAULIPS20100. DAT; PR=C£.PAULIPS2010013.LIS/N 
$ @L.PAULIRUN 

z Aaa, TRU »PAULIPS20100. DAT; PR=C.PAULIPS2ZO010014.L1IS/N 
$ GTSTRUDL C.PAULIPS20100. DATi PR=C.PAULIPSZ010015.LIS/N 
$ @L.PAULIRUN 

$ GTSTRUDL [.PAULIPS20100. DAT; PR=L£.PAULIPS2010016.LIS/N 
$ @CL.PAULIRUN 

$ GTSTRUDL [.PAULJIPS20100. DAT; PR=C.PAULIPS2010017.LIS/N 


oa wer oe woe oe wre wer wre ce ee we oe oe em cme wee Comm wee wet = wee Serr Cee oe cee Om C= come Cm ome we om ee oo om om ce oom Se ce cee Oe ee oe ee 2 ee we ae ee oe? ee eee oo ome ee eee ee ee ee ee eee eee ee oe 


a a cee Oe ee ee oe ee ee ee ee ee ee ee ee ee ws ws ee ee es ee es ee es ee ee es es ee es ee ee es es ee ee ee ee ee ee ee ee ee OE ae ee ee oe oe 


FOR C.PAVULIINC 
LINK INC 
RUN INC 
FOR iaecee OER 


ER 
[.PAULICREATE 
LINK CREATE 
RUN CREATE 
BEI ERC. +; + 
DEL CREATE. #; ¥ 
DEL LOADER. #i # 
DEL [£.PAULIPS20100. SAVi * 


tH Uh UH AF UF 2% th OF th OF th oF HF 
Tl 
S 
AN 





APPENDIX B 
SAMPLE GTSTRUDL OUTPUT 





GTICES 2.2  23-JSUL-1784 11:58:55.26 PASE 
pk 1.2 #8 Proprietary tc the Ceergia Tech Research Institute. 
erGTE STRUDL: GTIST8IO2. COB/IF_2=GIT_STRUDL: GTISTBIO2.OS/HODUMP/POOL_INCR=16334/1F 1=USERDAT .DS/ IF _7=fLOTFIL_DS 


DIRECT “L. PAULIPS201 00. SAY’ 


H-HEASS, Message number 08058057 
SREY=29.0, LSR=1.00, POSITION=45 A} 25.9 CTR’ 
| 


RRRAY RARER RRR RRRRK RRRWRK RRRKHK RRHRH RRNHA RENE RERKE RRRRRRKEEK 


# 

SH CITSTRUOL x 
HHH BCT GP 1783 - 83.92-759 # 
HH + 
HHRHH HMHHHH KKENH KH HK KEKHKK HH H 

HHASAHKHRY HHREHH HHHHHH HNKKRK RH HH KEXKHK Kx # 
HRUHNHRHHK HK HH HH HR HH NH HR HH HH + 
HERE HERE HH HHEKHH KH RH KK HM HK # 
MKNAH H HH ERE HH EKHNH NH KH HH HH HH # 
HEM RH HH HH HE KN RM EH HH HH HK ¥ 
HH HEM KHH HH KH XK RRHHHH HHHHH HHHHKK x 

#H HEH HH xy HH RH KENH  HHHHH KHRHHH x 
HH H 

HR QHNED BY ANB PROPRIETARY TO THE H 

HH CEQRETA TECH RESEARCH INSTITUTE 

H 


RRRAK RARKR RARE RRA RA RRA RRHRR RRNA KRHA RRR RR RRR RARRE REKRE A 


Tue WHITS - LERGTH hETCHT ANGLE TERPERATURE TIRE : 
MED TH BE THCH POUND §=RaDIGX = FAXRENHEIT © SECHND 
TE BEADING 

FENTARY TO FOOT HIGHWAY BRIDGE AvALYSIS PROJECT. 

METER STUDY IS USED TD INVESTISATE THE EFFECT OF 

HELE OX BRIDGE DECKS OF VARIDUS ASPECT RATIOS. 

"Ss. IN, DEG 

ERATE OFF 

DRDIRATES 

Y JOINTS ID 1,1 K 0.0 -16.41596683 ¥ 0.9 45.1052408 

TIHES ID ? % 48.9 

UPPORT 1 10 73 BY 9, ? TH 8i BY 9 

IRCIDERCES 

8 ELEMENTS ID 1,1 FROM 1,1 TD 19.1 1011.1 TD 2,1 

THES ID 8 FROW ? THY TO? TO? 

LEASES 

BY 2 HOMENT X BEAENT ¥ 

EY 9 MOMENT X, NORENT ¥ 

PROPERTIES 

TYPE ‘BPP’ THICKHESS 18.9 

S 

ALL 

9.15 ALL 

1 ‘FOOT Sil-y ROLLS ACROSS DECK CEXTER’ 
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2-0 43 





e700 
= -7 88 
2-77.59 
e-0.43 
Z 3.99 
e-1.99 
Z -0.98 
2 -2.10 
Z -2.88 
2 -2.82 
B-0.97 
2-90.83 
e -1.36 
2 -3.67 
2-27.99 
2-07.00 
Z -4.8/ 
2-90.53 
2-333 
y ANALYSIS 


) INFORMATION BEFORE RERUABERTAC. 


WA BAXDHIDTH IS 10 AND OCCURS AT JDINT 11 
SE BARDWIDTH IS 3.88? 
MAD DEVIATION OF THE BAXGHIDTH IS 2.846 


| REDUCTION HAS FAILED TO PRODUCE A BETTER HUAPERING. 


MUNBERING WILL PE USED. 


CONSISTENCY CHECKS FOR 64 RERSERS 6.38 SECONDS 


BARDWIDTH REDUCTION 3.99 SECONDS 
ENERATE «©6064 ELEAENT STIF. MATRICES = 12.54 SECONDS 
SSEMBLE THE STIFFRESS MATRIX 12.91 SECONDS 
RUCESS «= 81. «JINTS 3.66 SECONDS 
OLVE WITH 821 PARTITIONS 21.24 SECONDS 
ROCESS «=—«-«-B1s« SINT DISPLACERERTS 1.98 SECONDS 
ROCESS 64 ELEMERT REACTIONS 9.11 SECONDS 
ROCESS «©—64 ELERERT STRESSES 9.14 SECONDS 
STATICS CHECK 2.46 SECONDS 


AVERAGE PRINCIPAL BENDING RESULTANTS 
WE SPECIFICATION WISSING - HIDDLE SURFACE ASSURED 
HT LIST MISSING - ALL ASSUMED 


CTICES 1.2 
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PAGE 





HHRRR RRRHH RRRHH RRR RHE 


LTS OF LATEST ANALYSESS 
SHH MMH MM MC RE EH 


LEM - SREN=20. TITLE - HONE GIVEN 


VE UNITS IKCH KIP DEG DESF SEC 


m- 1 


CTICES 1.2 


ee ans ee ee om oe at aw ae em et ee ee ee et et ee es es ee ee et ee ee ee em ee ee) Se en en ne ee ee es ee ee ee ee we ee ee ee ee ee ee ee es et es ee es 
age amp ane om GD me mm me mm om om om et SD 00 wt et at nt nt et ns Dm m= at 8 et ee 8 em am me me me me mn me ne Se a es ee es ee ee ee a es ee De ee ee ot oe oe 


FDUT SU-4 RELLS ACROSS DECK CENTER 
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13-JUL-17 84 12:91:26.97 


rage 


ne ED ee Se ee? ee es ee ee ee ee oe ee ee ee ee et et ee ee) ee ee ee et et es oe et en Oe ee me ee ee OO Se Te ae Gan besa ee TS wee Or Ee SF 2" et ts ee ee ewe at ot Ge oa oF ee at ee eee ee eee a 


IPAL FENDIHG RESULTAXTS BASED ON AVERAGE BENDING RESULTANTS 


KUMEER OF ELEREXTS 


USED IN AVERAGING 


1 


2 


fd 


0. JF4OLIE+O9 


0. 2HO002E+09 


§.462039E+09 


0 GAPPBLE+99 


0. JOPTZSE+0O 


0. 368600E +09 


0. 619278E+90 


0. fASS0FE 109 


0.120179 E +01 


0 .116254E+91 


0. S020S7E+99 


9.377527E+09 


9.226925 +00 


0.126985E+09 


9. 182419E+00 


H2 


~0 7937 TEET)9 


“9 6463 249E101 


—0 8511632101 


~0 8707875401 


0 2353 20E191 


“Q .SSHOTPETOL 


9 FIO 70E+01 


~9 187 S4E+101 


9 HIL6H2E-02 


-9 U229USEtO1 


—) S867 J0EtOL 


-0 E667 90E401 


-) .PILALOE+O1 


-9 .825F16E101 


—O . S4SbW2EtO1 


HX Vii 


9.36 4875E 409 


§.3390626 +01 


0 WuA/S3E+0l 


0. WO7A7SE+O1 


0 RO IRSPEHOL 


0. SLOWSPE+O1 


0. 202U9PE 101 


O.117142E +01 


0. 62805HE +09 


§ LIPSOE +L 


9. F18573E +01 


0 .4S297SE +01 


0.48 7002E +01 


QO. HIPO07E +01 


0. FI1BN2E +01 


ret 


Q.1P7971E+02 
0. PASSO9E+O1 
). 7EROOUEHOL 
6. 62049 7E+01 
§.48lor 2Et01 
0. 21398 9E+01 
-). 280933E+01 
—O. 14139 1£ +02 
-0. 28579KE+02 
0.45902 SE 102 
0.2273 2E+92 
0.83229 SE+O1 
9. 39724CE +02 
0.42751 5E+01 


—0. 71616 5E+00 





0. 283081E+00 


0. UPSSSLE+00 


0.9 74702E +00 


0.1 58749£+01 


~9. APSHS6E+O0 


9. FO887SE+N0 


9. 92U320£+00 


-9 .FINSZHE+OO 


9. 608510E+00 


-9. SLIS8SE-01 


0. S88H1LE+00 


0. 150138£+01 


0. 164314E+00 


-0. 2798SIEt0L 


9. 429925E4+91 


-9.H1OZB1Et01 


-9.29T389E401 


-0. 249978E+01 


~9. 12639KE+O1 


0. 174997E+00 


0.162551 E+01 


0. 13560SE+01 


9. 201679£401 


0. F/1TSHGE+O1 


“9. U228S4E +01 


~0.476627E+01 


“0. 2964351 E491 


-9 .LPBOSSE+O1 


0. Y20514E +90 


“0 4276475401 


9. 228644E+01 


~9 L22598E +1 


“9 .171185E+91 


“9 .620562£+)1 


“9 F4ISZIE+91 


—0 1009 62E +02 


“0 .B9S776Et)L 


~9 7487 62E+91 


~) 9226232401 


—Q. 2984 SHEtOL 


~9 .135872E+)1 


~Q 2485 67E+09 


~9. 674211 +01 


-9 1995 95E192 


-0 .11é6885E+02 


“9 FSH2B9E +01 


“0.87089 7E+O1 


9 6184O7EtO1 


0 T4883 68E+01 


“9 .161967E+01 


-Q.125626E+01 


~9. 616671£+01 


“9 10059 7ET92 


—¢ 1189 FHE+02 


~9 .193074E+02 


~9 .86479GE+01 


-0 .670112E+01 


“0 .FS61PSE+01 


GTICES 1.2 


0.228227E +01 


0.137105E +04 


0. 1O0P789E+01 


0. 16P9S5E +01 


0. 270008E +01 


0.426368 +01 


0. 458592E +0} 


0. NO 2S59E +01 


0. FUAPSSE +01 


0. 2Z57254E +01 


0.178638E +01 


0. LY FOOSE +O 


9. 20S440E +00 


9. U9 8SHhE HL 


0. 3180636 401 


0.37 9289E +01 


0. F40S4/E +01 


0. ARFU6OE OL 


0. 2U5706F +01 


0. 183184E +01 


O.1S1P60E +01 


0. 129516E +91 


0.21 7U96E+O1 


0. 34221KE +01 


§.333236E+01 


0.42 7057E +01 


0. 284220£ +01 


0. ZI6027E +91 


0.15 7067E +01 


-0, 89050 6E+01 
~9 22748 2E402 
~} 7 L271 SE 402 
9 4G2509E+02 
). 15542 £402 
). PFEZI9E+OL 
>. SHudSBEtOL 
0. 47S80hE+00 
~}. 28729 5E+01 
~9 11060 2E+02 
~O. 2UPIZHE+O2 
—? QUPS2AE+02 
0.526 79E402 
0. $8692 5E+01 
0. 4572Z25E +01 
0. 22803 HE FOO 
9. TIO200E+01 
-9 .DVSVLZE+O1 
-0. 1F3459E+02 
0.27021 5E+02 
0. 452945E+02 
-O USI2UPE+02 
9. 2039S0Et01L 
0. 116965E+01 
O. 2u241GE+O1 
-0. 69301 FE #00 
-). 38892 6E+O1 
-0. S150 9E+O1 


0. 24445 7Et02Z 
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PA 


6 


E 








0. Lh0617E+91 


0. 208701 E01 


~9.120304E+91 


~) F6ITGE+O1 


“9 USTIPOE +01 


~0 U29291 E401 


9 A256PFE+OL 


~0.270676E+01 


~O .LU2403E+01 


0.848907 £409 


0. 212H1KE+O1 


0. 4O176SE 4990 


“9. 687197 E400 


~ 29 ,142582F +01 


~0.1742935E+01 


~9. 1642452401 


~0.119085E+01 


“0.7 71ISP7E+00 


0.225274 E+99 


9.125957E+91 


0. 994257E+00 


0. 217425E499 


0. J51920E-01 


~0.127372E409 


~0.207859E+09 


—9, LOUPTLE+O9 


~9.1L4788hE-91 


9. 2bu7FE+O0 


0. LIGSPOE+OL 


“0 .LUUP PPE +01 


~0 1985 7BE+O1 


-) 5027 67E+01 


~9 .191225E+02 


-0 .117610E402 


“O LOPOLHETO2 


“0 BYP245EtO1 


“0.72501 7Et01 


oa 


. 438322E+91 


.PY18 BLE+00 


t 
© 


~0.187712E+)1 
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APPENDIX C 
CONTOUR PLOTS 
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FIGURE C.1 CONTOUR: MAJOR PRINCIPAL 
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FIGURE C.5 CONTOUR: Z DISPLACEMENT FOR CENTRAL 10 KIP LOAD, SKEW = 0° 
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FIGURE C.7 CONTOUR: Z DISPLACEMENT FOR PEAK SU-4 CENTER LOAD, SKEW = 0° 
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